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Sommaire. — Quand le champ H croît, l’aimantation # d’un ferromagnétique cubique polycristallin 
tend vers Ta saturation Y, selon la loi d'approche ; 


2) a) 


et une théorie actuellement admise montre que b est égal à 8 K2: 105 3? où K est la constante d’aniso- 
tropie. L'auteur montre que cette théorie n’est pas correcte car elle a été obtenue en négligeant les inter- 
actions magnétiques entre les cristallites dont le rôle est considérable. Les vérifications expérimentales 
qu’on en a données sont tout à fait illusoires, soit par suite d’interprétations trop sommaires, soit par 
l'emploi de substances à tensions internes résiduelles trop grandes. 

L'auteur développe ensuite une méthode nouvelle de calcul qui permet de tenir compte des inter- 
actions magnétiques et obtient ainsi une loi d'approche rigoureuse qui n’a pas la forme (A). Cependant, 
dans la région pratiquement la plus intéressante, c'est-à-dire de 200 à 1000 œrsteds, elle peut se mettre 
approximativement sousla forme (A), mais la valeur de b est à peu près deux fois plus petite que la valeur 


donnée plus haut. De nouvelles mesures effectuées sur du nickel bien pur et bien recuit ont confirmé ces 


résultats. Les formules obtenues permettent de déterminer X au moyen de mesures effectuées sur des 
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échantillons polycristallins. 


te Introduction. 


1. Détermination de la constante K d’aniso- 


tropie, par l'étude de l'approche à la saturation. 


— Cas d’un monocristal. — On sait que l’énergie 
libre F d’un cristal ferromagnétique est la somme de 
plusieurs termes dont l’un F,, appelé énergie magnéto- 
cristalline, dépend de l'orientation de l’aiman- 


tation spontanée Y, par rapport aux axes cristallins. 


- Dans un cristal cubique, l'expression suivante 


donne une représentation satisfaisante de F,,, en 


fonction des cosinus directeurs «, B, y, de #, par 
rapport aux axes quaternaires 


Pre (a p2+ B2y2 + y202). ( 1) 


(:) Les travaux expérimentaux et théoriques qui font 
l’objet de cet arti le ont été poursuivis et achevés pendant 


_ lPoccupation. La parution en a été retardée par suite de 


difficultés matérielles. Devant les retards d'impression qui 
s’annonçaient déjà en 1945, une brève note a été publiée 
sur le même sujet (C. R. Acad. Sc., 1945, 220, p. 814). 


Si le cristal est soumis à.un champ uniforme Hi, 

l’énergie libre comprend en outre un second terme F, 
> 
Dali Ha. 

Dans un cristal parfait, la position d’équilibre 
de ë s'obtient en écrivant que l'énergie 
libre F—F,<+F; est minimum. Lorsque le 
champ H devient de plus en plus intense, tout en 
gardant une direction donnée par rapport au 
cristal, l’aimantation spontanée tend à devenir 
parallèle au champ et sa composante # suivant 
cette direction, tend vers l’aimantation à satu- 
ration 3. On montre facilement que l’approche 
à la saturation s'effectue suivant une loi de la forme 


T = Y G KE 
d' = RÉÉTD , ( 


où le coefficient b dépend de l’orientation de H 
par rapport aux axes du cristal. 


RO 
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Cas d'un polycristal. — La formule précédente, 
applicable à un monocristal, n’a fait l’objet d'aucune 
vérification expérimentale, mais depuis longtemps on 
s'est proposé de l'appliquer à un polycristal. Le 
premier, Akulov (2?) a calculé la valeur de b relative 
à un polycristal en faisant la moyenne des valeurs 
obtenues en donnant à un monocristal toutes les 
orientations possibles; on obtient ainsi 


Be (3) 


Gans () et plus tard Becker (?) ont repris le même 
raisonnement et obtenu le même résultat. 


Vérifications expérimentales. — Entre temps, 
Czerlinski (5) a montré expérimentalement que les 
polycristaux de fer et de nickel obéisssaient à une 
loi d'approche de La forme (2) et, à partir des valeurs 
de b ainsi obtenues et compte tenu de la formule (3), 
il en a déduit des valeurs de X qui, entre — 150° 
et 1509 C, concordent avec les valeurs de la litté- 
 rature déterminées par les méthodes habituelles 
applicables aux monocristaux. De son côté, Polley (5) 
a repris l’étude du polycristal de nickel et montré 
qu’en réalité il fallait remplacer la loi d'approche (2) 
par une doi plus complexe de la forme 


7 b 
= T4 — > — —— Y IT =: À 
d (i 7 7. + x 1H; (4) 


. °«d ù 
mais, en supposant que le terme en H provenait 
d’un mécanisme inconnu mais n’altérant pas le 

DR | : 
terme 2» il a pu, comme plus haut, calculer à 
différentes températures des valeurs de’ K qui 
s’accordent encore avec celles de la littérature. 

Il semble donc que nous soyons en possession 
d'une méthode éprouvée de détermination des 
coefficients d’anisotropie K, à partir de mesures 
efféctuées sur des polycristaux, ce qui présente un 
intérêt expérimental évident. Malheureusement,’ on 
peut élever contre cet ensemble de résultats de 
graves objections, tant du point de vue théorique 
que du point de vue expérimental. 


2. Objections d'ordre théorique à la méthode 
précédente. — Il ne semble pas légitime de calculer 
le coefficient b relatif à un polycristal en prenant 
la moyenne des valeurs de b relatives aux cristallites 
élémentaires supposés isolés, car on néglige ce 
faisant les champs démagnétisants internes déve- 
loppés par les charges magnétiques qui proviennent 
de la non-uniformité de l’aimantation du poly- 
cristal. 


OAKRULOV 227 RPRUISIL 1000 0060 SD MS D: 

(5) GANSs, Ann. der Physik, 1932, 15, p. 28. 

(1) Becker et D6RING, Ferromagnetismus, p. 168 (Springer, 
Berlin, 1939). 

(5) CzerriNsKkt, Ann. der Physik, 1932, 13, p. 80. 

(f) Pozrey, Ann. der Physik, 1939, 36, p. 625. 


ES 


Pour un cristallite isolé, la présence a Mb 


démagnétisant n’introduit pas de difficulté spéciale : : 


le cristallite est en moyenne sphérique de sorte que 


son coefficient moyen de champ démagnétisant est “ 


indépendant de la direction et égal à . D'autre 


3 
part, comme la théorie suppose implicitement un 


champ appliqué suffisamment intense pour qu'il 
n’existe qu’un seul domaine élémentaire par cris- 


tallite, il en résulte que l’aimantation d’un cristallite 
est toujours égale à l’aimantation spontanée et que 


les seuls changements d’aimantation à considérer. 
proviennent de la pure rotation d’un vecteur de 


grandeur 4, invariable (7, rotation au cours de 
laquelle l’énergie de champ démagnétisant reste 


I 
invariable et égale à — 1% 42 FerOS CHE 


DS: 


donc pas à tenir compte de ce champ démagnétisant. 


Il n’en est plus de même pour l’ensemble des 
cristallites juxtaposés qui constitue le polycristal. 
Supposons en effet que les cristallites possèdent 
ainsi réunis, le même état d’aimantation qu’à l’état 


isolé. Comme la direction de l’aimantation spontanée 


varie d’un cristallite à l’autre, il apparaît sur les 
surfaces de séparation des charges magnétiques 
fictives qui créent des champs magnétiques et 
augmentent ainsi l’énergie libre. Mais, tandis que 
dans un cristailite isolé cette énergie conserve la 
même valeur lorsque l’aimantation devient parallèle 
au champ appliqué, ici au contraire cette énergie 
supplémentaire disparaît à la saturation lorsque 
tous les cristallites prennent des aimantations 
parallèles (on suppose que le polycristal a la forme 
d’un barreau allongé dans la direction du champ). 


Il en résulte nécessairement des forces qui tendent 


à faciliter la saturation. 

On peut apprécier l'effet de ces interactions par le 
raisonnement suivant : Si l’on extrait d’un poly- 
cristal un cristallite élémentaire, le champ qui 
règne dans la cavité correspondante est, en moyenne 

fé 
égal à H + 373, en désignant par à l'aimantation 


moyenne du polycristal. C’est ce champ total qui 
agit eflectivement sur le cristallite. Comme au 
voisinage de la saturation, # est pratiquement 


égal à %,, il faut donc remplacer le terme en 
b 
4 

(a + = 3 Tr J ) 

+ est voisin de 500 C. G.S., 


VE 


de la loi (2) par un terme en + Dans le 


cas du nickel, où # 


et dans un champ de 1000 œrsteds, on constate que - 


le nouveau terme est neuf fois plus petit que le 
terme primitif. Sans doute le raisonnement esquissé 
ici est-il trop grossier pour nous donner un résultat 


(*) Pour cette question, se reporter au mémoire de NÉEL, 


Les lois de l’aimantation et de la subdivision en domaines . ” 


eo ce d’un monocristal de fer, J. de Physique, 1944, 
5, p: 241-2517. 


Il n'y a 


It 


pq 
y 


l 


Pos mais il suffit à montrer que les interactions 


entre cristallites sont considérables et qu’il ne peut 
être question de les négliger. F 

. Finalement la loi (2) avec la valeur de b donnée 
en (3) est grossièrement inexacte. 


. 3. Objections d'ordre expérimental. -— Il 
_ s’agit alors d’expliquer dans quelles conditions les 
_ auteurs précités ont pu obtenir des valeurs correctes 
de X à partir de formules inexactes. 


Expériences de Czerlinski. — En ce qui concerne 


_ ces expériences, une discussion serrée du mémoire 
_ de l’auteur, qu'il n’y a pas intérêt à reproduire ici, 
. montre que ses résultats sont beaucoup mieux repré- 
_ sentés par une formule du type 


Bt en) “à 
ls s ter 2 la f ] ‘il indique 
ans terme en ;%» que par la iormule qui indique, 


« . 
ne comprenant pas de terme en ;;; mais avec un 


à 


b : DRE 
» et qui présente des déviations systéma- 


terme en 
FREE 


tiques par rapport à l'expérience. On obtient d’ailleurs 


selon cette nouvelle interprétation une valeur de a 
égale à 5,2 qui est du même ordre de grandeur que 
les valeurs obtenues par d’autres auteurs pour du 
fer pur (%). Aiïnsi les expériences de Czerlinski, 
loin d'appuyer la théorie d’Akulov-Gans-Becker, 
montrent plutôt que la valeur du terme en _ 
prévue par celle-ci est trop grande. Il faut cependant 
remarquer que la présence d’un terme supplémen- 


taire en 7 dû probablement aux imperfections 


des cristallites (5), rend l’interprétation incertaine 
car il est à craindre que ces imperfections ne se 
manifestent aussi par un terme supplémentaire 
AO: > | 

.en T° 

Expériences de Polley. — De ce point de vue, les 
expériences de Polley sur le nickel se présentent 


plus favorablement, car le terme en F y est très 


faible, mais l'échantillon de cet auteur doit néanmoins 
être encore assez éloigné de la perfection, comme 
le montre la méthode suivante. 


4. Signification et calcul du travail d'aiman- 
tation. — Dans l’état actuel de nos connaissances, 


nous ne savons pas calculer la courbe d’aimantation. 


d’un polycristal, mais il existe une grandeur dont 
la détermination expérimentale est facile et qui, 


(5) Néez, La loi d'approche en 1 : H, J. de Physique, 1948, 
9, p. 184. 
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indépendamment de toute théorie, peut être calculée 
en fonction de la constante d’anisotropie : il s’agit 


du travail d’aimantation. 


Du point de vue expérimental, le travail d’aiman- 
tation est égal à l'aire comprise entre la courbe 
d’aimantation (aimantation en fonction du champ), 
l'asymptote à la saturation et l’axe des ordonnées. 
Cette méthode suppose qu’il n’y a pas d’hystérésis: 
on peut cependant encore, comme l’a montré 
Kersten (?) dans un cas analogue, évaluer le travail 
réversible d’aimantation quand l'hystérésis est 


faible : en mesurant l’aire comprise entre la branche 


descendante du cycle limite et l’asymptote à la 
saturation (fig. 1) car cette branche SR du cycle 
est sensiblement réversible. 


O AE 


Fig. 1. — La branche descendante SR du cycle limite est 
sensiblement réversible : l’aire hachurée correspond donc 
au travail d’aimantation. 


Du point de vue théorique, le travail d’aimantation 

est égal à la différence entre l’énergie interne à la 
saturation et l’énergie interne dans l’état initial, 
quand l’aimantation est nulle. 


Dans cet état initial, l’aimantation spontanée 
des différents domaines élémentaires est dirigée 
suivant l’une des directions dites privilégiées, corres- 
pondant à un minimum de l’énergie magnétocris- 
talline. Lorsque X est positif (cas du fer), le minimum 
de l’énergie magnétocristalline F,, est zéro; lorsque K 
est négatif (cas du nickel) ce minimum est égal 


UK , ; : É À 
à - + D'autre part, à la saturation, l’aimantation 
re) 


spontanée est parallèle au champ appliqué, de sorte 
que la valeur moyenne de l’énergie d’un polycristal 
est égale à la valeur moyenne de F, donnée par 
l'expression (1), quand le vecteur unitaire «, $, y 
prend. toutes les orientations possibles : un calcul 


E 


élémentaire montre que cette moyenne est égale à LR * 


Finalement, il en résulte que le travail d’aimanta- 
tion d’un polycristal est égal à : ergs : cm°, quand K 
K — 22 


Re ne LIKE > K ; : 
est positif, et égal à À Ur cg: CM; 


quand À est négatif. 


() KERSTEN, Z. Physik, 1932, 76, p. 505. 
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Application aux expériences de Polley. — Or, 
d'après les données de la littérature relatives aux 
monocristaux de nickel (cf, Becker, Loc. cit, p. 123), 
K est voisin de —5,0.r104ergs : cm*, d’où il résulte 
qu'à la température ordinaire, le travail d’aiman- 
tation d’un polycristal de nickel est voisin 
de 6,6.10% ergs : cm. 

D'autre part, d’après les données numériques et les 
graphiques du mémoire de Polley, il résulte qu’à 14° C 
le travail d’aimantation du nickel qu'il a étudié 
était voisin de 13.10% ergs : cm”, valeur double de la 
valeur théorique. Malgré le recuit qu'il avait subi, 
ce nickel était donc probablement affecté de tensions 
internes relativement considérables altérant de 
façon importante la courbe q’aimantation; on sait 
en effet que, par suite de la grandeur de sa magnéto- 
striction, le nickel est particulièrement sensible 
aux tensions internes. 


5. Nouvelles mesures sur le nickel. —— Dans ces 
conditions, il était intéressant de reprendre la mesure 
du coefficient b sur du nickel pur, débarrassé de 
tensions internes. 


Méthode. — À cet effet nous avons étudié, à la 
température ordinaire, un fil de nickel pur. Le fil, 
de o,, mm de diamètre et de 130 cm de longueur, 
était tendu à l’intérieur d’un solénoïde de 95cm 
de longueur. La bobine de flux, convenablement 
compensée de manière à indiquer directement au 
galvanomètre balistique les variations de 47, 
avait 20 cm de longueur. Toutes dispositions étaient 
prises pour pouvoir tracer des cycles, point par 
point, à partir d’une origine arbitrairement choisie. 
Le fil avait été recuit 3 h à ro000 €, dans de l’hy- 
drogène soigneusement purifié. Des recuits ultérieurs, 
poussés jusqu’à 1200° C, n’ont pas modifié de manière 
sensible la forme de la courbe d’aimantation. 


Résultats. — Ils sont donnés par les figures 2 
et 3. Le champ coercitif de ce nickel était très faible : 
égal à 0,66 œrsted. On pouvait donc espérer qu'il 
était bien libre de tensions internes. Nous avons 
en effet trouvé que l’aire comprise entre la branche 
descendante du cycle limite et son asymptote était 
égale à 6,4. 10% ergs : cm”, à la température de 279 C; 
cétte valeur est extrêmement voisine de la valeur 
théorique obtenue plus haut. 

Pour préciser la loi d'approche à la saturation 
du nickel ainsi traité, nous avons sur la figure 3 
représenté l’aimantation en fonction de ES Entre 150 
et 750 œrsteds, on obtient ainsi sensiblement une 
droite qui s’extrapole vers une aimantation à satu- 
ration égale à 497,5. Dans cet intervalle de champ, 


les termes en À et en ;H du développement (4) 


ne sont pas sensibles. La formule (2) suffit alors 
pour représenter les expériences avec un coefficient b 


. 


“N° 6: 


égal à 290, à la température de 270 C: Pour cette 
même température, en interpolant les résultats 
de Polley, on obtiendrait b — 550, valeur sensi- 
blement double de la précédente. 
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Fig. 2. — Courbes  d’aimantation dans des champs moyens 
de différents échantillons de nickel : A, fil tel qu’il sort 
de la bobine (H, = 8,7 œ); B, le même fil recuit 15m 
à 9oo° dans l'air (H, — 5,4 œ); C, le même recuit 3h 
à 1000° dans l'hydrogène (AH, = 0,7 œ); D, fil étudié par 
Polley. 
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Fig. 3. — Approche à la saturation dans les champs forts, 
en fonction de l'inverse du carré du champ : C, fil de nickel 
rvecuit 3 h à 1000° dans l'hydrogène; D, fil étudié par Polley 
(d’après ses mesures). 


Cette différence entre les propriétés de deux échan- 
tillons de nickel pur ayant subi des traitements 
thermiques voisins est difficiie à expliquer : remar- 
quons simplement que l'échantillon de Polley 
était plus massif que le nôtre, ayant 5 mm de 


_ diamètre, et que les impuretés pouvaient y avoir 
plus de difficultés pour diffuser à la surface et 
s’éliminer; en outre, Polley ne semble pas avoir 
À pris de précautions spéciales pour purifier l’hydro- 
gène. 
Quoi qu’il en soit, il en résulte que du point de 
. vue expérimental, pour du nickel aussi bien débar- 
_rassé que possible des tensions parasites, le terme 
en . de 
et 700 œrsteds, deux fois plus petit que le terme 


la loi d’approche est, entre 200 


théorique, calculé en négligeant les interactions 


entre les cristallites. 


II. Théorie de la loi d'approche 
d'un polycristal. 


6. Représentation analytique des forces 
magnétocristallines dans un polycristal. — 
Exposons maintenant une méthode de calcul appli- 
cable à un polycristal constitué de cristallites élémen- 
_ taires, dont les orientations sont réparties au hasard, 
placé dans un champ magnétique suffisant pour 
qu’en tout point l’aimantation soit voisine de la 
saturation. En chaque point, l’aimantation est en 
équilibre sous l’action du champ magnétique 
appliqué Æ et des forces magnétocristallines. Tant 
qu’on se borne-à étudier l’approche à la saturation, 
ces dernières forces sont équivalentes du point de 
_ vue formel à un champ magnétique F, d'intensité 
convenable, perpendiculaire au champ H : naturel- 
lement la grandeur et l’orientation de F changent 
d’un cristallite à l’autre. 

Découpons d’autre part dans le polycristal un 
_ cube de côté L, très grand vis-à-vis des dimensions 
des cristallites. Il est clair qu’en juxtaposant les 
uns aux autres, par des translations suivant trois 
directions rectangulaires, des cubes identiques au 
premier, on obtient un polycristal périodique indéfini 
de propriétés pratiquement identiques à celles du 
polycristal initial. Ceci posé, rapportons ce poly- 
- cristal périodique à trois axes rectangulaires OXVZ. 
La composante du champ F suivant une direction À 
quelconque est développable en une série de Fourier 


= Ÿ Myurexpi 2 (pe + aq +rs)) (6) 


où le signe somme est étendu, ici comme dans tout 
ce qui va suivre, à toutes les valeurs entières de pgr 
de —c à + oc. Les coefficients M, sont des 
quantités complexes qui, F; étant réel, satisfont 


à la relation à 
ME pq, Ga Mr ( £ ) 


en désignant par un astérisque * l'imaginaire 
conjuguée. 

En découpant le cube initial dans un autre endroit 
du polycristal, on obtient naturellement des coef- 
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ficients M,,,. entièrement différents de sorte qu'il 
importe plutôt de connaître leurs valeurs probables 
que leurs valeurs réelles pour un développement 
particulier. Comme le polycristal est en moyenne 
isotrope, il en résulte que, comme nous l’avons 
déjà montré sur un exemple analogue (#), la loi 
de probabilité qui fixe la valeur de M,,, ne dépend 
que du seul paramètre : o—p?+ q?+r?, car le 
développement doit conserver la même forme pour 
une rotation quelconque des axes. En particulier, 


- la valeur moyenne M,,,, M*,,, du produit M,,: M*,,r 


est une fonction positive de o. 

Supposons alors les axes orientés de manière à ce 
que OZ soit parallèle au champ appliqué H; les 
composantes F, et F, de F suivant OX et OV 


sont représentables par deux développements du 


type (6), identiques en moyenne par raison de 
symétrie. Les coefficients correspondants, M,,; 


et M',,,. obéissent aux relations suivantes, où l’on 


a pris la valeur moyenne, par rapport à toutes les 
positions possibles du cube initial, de la quantité 
située au-dessous du trait 


M 0 


par 


Mogr My MopMour= En (8) 


où }, est une quantité positive qui ne dépend que 
de ©. En outre, puisque les M et les M’ sont indé- 
pendants, on a 


M pqr Mgr = 9; M3» Mpqr = 9: C9 
7. Loi d'approche sans interactions entre les 
cristallites. — Désignons par « et 6 les angles 
avec OZ des projections sur les plans XOZ et YOZ 
de l’aimantation #, en un point. Puisque les angles « 
et 8 sont petits, la conditions d’équilibre de l’aiman- 

tation s'écrit simplement 

F, 
B= 75 (10) 
et, au second ordre près, l’aimantation macro- 
scopique #, égale à la valeur moyenne de la projection 
de l’aimantation sur OZ, s’écrit 


a + (2 
ETS a , : =; Ÿ s\I 


Pour obtenir F}?, il faut, après avoir élevé l’expres- 
sion (6) au carré, effectuer successivement les 
moyennes par rapport à æyz et par rapport à toutes 
les valeurs possibles des coefficients M; on obtient 
finalément, compte tenu des relations (7) et (8) 


F2 F? 
Æ AE 11 
2 H? ) (ED) 


FF? = V My HS F,, (19) 
d'où la loi d'approche 
= 91 ne, (13) 
FR ) 
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Comme cette loi doit être identique à la loi 
donnée par les méthodes classiques, on déduit, 
après avoir rapproché les relations (2), (3) et (13) 


Sr, 8 K? 
DRE 


8. Introduction des interactions. -- Nous 
allons maintenant introduire les interactions entre 
les cristallites qui avaient été négligées par Akulov 
et ses successeurs : Gans et Becker. Ces interactions 
. se manifestent par l’existence d’un champ démagné- 
tisant interne dérivant d’un potentiel V et provenant 
de la non-uniformité de l’aimantation du poly- 
cristal. La condition d’équilibre (10) doit être 
remplacée par la suivante : 


1 ov D 0: 
rte eh Perte) 


D'autre part, le potentiel V satisfait à l'équation 
de Poisson : AV + 4{rm—o avec m —= — div à. 
Mais comme les composantes locales de l’aimantation 
sur les trois axes sont, au second ordre près, respec- 
tivement égale à « ,, 8 #, et 3, on en déduit 
Éd TO 


(14) 


(15) 


si 1 dF GER A 2 16 
Ho ne Le) 
et l’équation de Poisson prend la forme 
(r+ ce) cu DV a 
A }\Vor ae] "07 
:Ard,/dF. or. ds 
Lie ai . 


Des développements du type (6) qui définissent F> 
et F,, on déduit 


or oF 217 
x DC M pqr GE qM Dyr ) 


2 
à (18) 


ÿ Rd 2 © 
+= — F e? : SARG 
ir . AU au foparar= TE Je Fe des 


où G est l’intégrale définie 


T 


cf SE Re d 
À (+ b sin?o} Aa 


qu’un calcul élémentaire donne égale à 


(24) 


I 
FR 2(b +1) 


PR 2 are RENE 
2(b +1 B: CAB CE RE MES 


D'autre part, la même transformation de somme 
en intégrale appliquée à la somme È2F, donne 


co 
Die x f ef, de. 
0 


(26) 
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L’équation (17) est donc satisfaite par un un dévelop- 2 | 
pement de la forme ” | 


pds EGpe + gr + rs) (19) 4 


où les coefficients w sont donnés par 


bL(pMpqr+ 4 Mur) 
2ir[O(p+g)+p+g+r] 
EDR 


Ppqr = à D | 
(20) à) 


avec _ D =: 


En reportant cette expression de V dans la . 
relation (15), on en déduit comme valeur moyenne … 
de x? par rapport à xyz 


PR Ep(pM+qM) | 
TP di ee gi}? 
nr" | 


bp(pM+qM*) 
b(P?+ g°)+ 6° 
et une expression symétrique pour f2. 
En effectuant ensuite la moyenne par rapport à « 
toutes les valeurs possibles des M et des M’ et « 
compte tenu des relations (8) € et (9), on obtient 
comme valeur moyenne & + B? de ao? + f2 


ep Mr JS 


Comme le nombre des termes compris sous 
le signe somme est très grand et que la proba- 
bilité d'en rencontrer un dans le volume dp dq dr 
est partout simplement égale à dpdqdr, on … 
peut remplacer le signe somme par une inté- = 
grale et en effectuant le changement de variable : 
p=?psinocost, g—psinvsint, r—pcos9, 4 
dp dg dr = p?sin © do do dû, on constate que les « 
variables se séparent et en intégrant par rapport . 
à 0, il vient 


œ* 


1) 


(22) 


(23) 


Finalement en substituant l'expression (23) dans É. 
l’équation (11) et compte tenu des relations (26) 
et (14), on aboutit à la loi d'approche suivante 
où Gest donné par l'équation (25) avec b égal . 
à 4TYS 4 
TR 
(27) 


sfr 8 k? A 


1058712 9 


9. Discussion de la loi d'approche (27). — 


Supposons d’abord que le champ appliqué soit très. | 


grand; dans ce cas b est très petit et l’on trouve qu’en 


négligeant les termes du second ordre en b, _ se 


,* La loi (27) devient alors 


se 
ak ) 
5 g 2 
RPM # . (28) 


105 42 (u Le = so) 


: On constate alors que les interactions entre les 
cristallites ont pour effet d'ajouter au champ appliqué 


un champ égal à À 7 Y,; on reconnaît là le terme 


. classique de Lorentz qu’un raisonnement sommaire 

nous avait donné plus haut. Mais, du point de vue 
_ pratique, il est impossible de se placer dans des 
. conditions où la formule (28) soit valable : en effet, 

comme pour le fer 473, est déjà égal à 21 500, 
il faut atteindre 100 000 œrsteds pour que b, à la 
rigueur, puisse être considéré comme petit. On 
_ sait au contraire que, ee l’étude expérimentale 


de la loi à ‘approche en ,;; les champs utilisés re 


_ 
dépassent pas 2000 œrsteds. Dans ces conditions, 
b doit plutôt être considéré comme grand. Le calcul 


_ numérique de l'expression (25) montre que G :», 


z. x I At à EE Er 
égal à 0,5 quand ; est nul, croît à peu près linéai- 


: SC à Mare è 
rement, avec 5 Jusqu'à ; — 0,3. C’est ainsi qu'on 


‘obtient les chiffres du tableau suivant : 


TABLEAU. 
I A 
pt‘ 0,02 0,050 0,116 0,270 
RU À É 
PRE 0,506 0,513 0,531 0,576 


en se limitant au domaine accessible à l'expérience. 
Dans cet intervalle, G varie assez peu avec H pour 
que la loi d'approche réelle puisse être confondue 
en is approximation avec une loi rigoureuse 
ne conformément à l'expérience, mais le coef- 
ficient du terme en _ est sensiblement deux fois 
plus petit que celui que l’on obtient en négligeant les 
interactions. 
‘Dans les mesures sur un fil polycristallin de nickeï 


d'étude de la loi d'approche en 
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qui ont été exposées précédemment ($ 5), le domaine 


I GR ‘ 

mn S étendait 
SPAS [ pote 

de -H= 150 4H = 550 œrsteds :: ; variait donc 

de à 102% à o,1 et la valeur moyenne correspondante 


tie & d’après le tableau, était d'environ 0,515. 


Le coefficient du terme en ayant été trouvé 


I 
12 


égal à 290, on en déduit d’après la formule (27), 


K = /,3.10tergs: em, pour {—9270C. A ‘cette 
même température, les mesures de Brukhatov 
et Kirensky (1) sur un monocristal donnent 
K = 4,16.104 ergs: emÿ. L'accord entre ces deux 
valeurs est excellent. 


10. Conclusions. — Ainsi, en ce qui concerne 
le nickel, l'introduction des interactions entre les 
cristallites conduit à une formule qui s'accorde 
avec les résultats expérimentaux. 

Inversement, il semble qu'on puisse utiliser a 
formule (27) pour déduire de la loi d'approche la 
constante d’anisotropie K; à titre de première 
approximation, on pourra poser G — 1, c’est-à-dire 
utiliser la formule d’Akulov-Becker, mais avec un 
coefficient deux fois plus petit. 

On peut également déduire K de la mesure de 
l'énergie d’aimantation, en suivant la méthode 
exposée dans le paragraphe 4, mais cette méthode 
paraît moins sûre et moins commode. 

D'une manière générale, il sera bon de s’assurer, 
en utilisant les rayons X par exemple, que la 
substance est réellement isotrope et que les axes 
cristallins sont effectivement répartis au hasard, 
ce qui n’est pas toujours le cas des fils étirés. Il 
faut, en outre, que la substance soit aussi bien 
débarrassée que possible de tensions internes et ne 
contienne ni cavités ni inclusions non magnétiques : 
en effet, comme nous l'avons montré antérieu- 
rement (8), la présence des cavités perturbe profon- 
dément la loi d'approche. Ces précautions étant 
prises, il semble qu’on puisse tirer un bon parti 
des mesures sur les polycristaux, dans l’étude de 
l’anisotropie magnétique des substances cubiques. 


(1°) BRuKHATOV et KIRENSKY, Physik Z. Sowjetunion, 1037, à 


12, °p-:002° 


Manuscrit reçu le 3 mars 1948. 
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SUR LE CALCUL DES BASES RAYONNANTES 


Par MAURICE PARODI et GEORGES PIRCHER:. 


Sommaire. — Dans cette étude, les auteurs donnent, en s'appuyant sur les propriétés du produit 
de compo$ition, d'usage courant en calcul symbolique, une méthode de calcul de bases rayonnantes 
ayant des maxima secondaires très faibles vis-à-vis du maximum principal. 


De nombreux travaux ont été effectués en 
Acoustique et en Radioélectricité en vue de déter- 
miner des bases d'émission ou de réception ayant, 
dans une direction de plan déterminée, un maximum 
principal aussi aigu que possible et des maxima 
secondaires négligeables devant ce dernier, parmi 
lesquels il faut citer ceux de Hs (2) et de Shel- 

kunoff (?). 

Nous nous proposons de ont que le calcul 
symbolique permet, en ayant recours à la notion 
de produit de composition, de donner une solution 
du problème. 


1. Introduction. — Considérons une base d’émis- 
sion linéaire résultant de la distribution continue 
sur une longueur finie ou infinie de centres d’émission 
vibrants en phase (fig. 1). 


Fig 1. x 


L’amplitude de la vibration résultante à l'infini 
dans une direction faisant un angle B avec la normale 
à la direction du groupement (ou directivité de la 
base) a pour expression. 


ARNO EUR ee 
fer rende) (2) 
0 


f(æ) étant l'amplitude de la vibration émise au 
point d’abscisse x (fonction de répartition), cette 
dernière étant comptée à partir de l’une des extré- 
mités de la base, À la longueur d’onde et j le symbole 
des imaginaires. 

Si l’on compare l'écriture (1) à celle qui donne 
le module de la transformée de Laplace + (p) de la 
fonction f (x) 


eue] f erpayaal 


| 
à 
# 
5 
e 


A(B) = 


() Akustische Zeils., 1941, 6, p. 9o. 
(*) Bell. Syst. tech. Journ., 1943, 22, p. 80, 


il apparaît que À (8) est égale au module de + (p) 
: PR 
quand on y fait p——)-snmé$. 


< 


Les tables de correspondances symboliques per- 
mettent done de déterminer la directivité d’une base 
linéaire ayant une fonction de répartition donnée. 


Nous pouvons remarquer que la fonction de répar- 


tition f (x) peut, au lieu de caractériser l’amplitude e 


au point d’abscisse æ d'une distribution linéaire, 
représenter la largeur comptée perpendiculairement 
à Ox, d’une base plane continue dont tous les points 
vibrent en phase avec la même amplitude; dans ces 


conditions. 
e(—: "= sin 4) 


représente la directivité de la base dans la direction 
de plans perpendiculaires à sa surface et de traces 
parallèles à Or. 


Soit, par exemple, une base plane rectangulaire 
de largeur b et de longueur a (fig. 2); f (x) est nulle 
quand x est extérieur à l'intervalle (0, a) et égale à b 
quand x est intérieur à cet intervalle. 


A(B}= 


Fig. 2 


Proposons-nous de déterminer la directivité A0 
de cette base dans le plan XO7Z. 


La fonction / (x) a pour transformée de Laplace 
L = er? « 


1? 


b 


et la directivité s'écrit, d’après ce qui précède, 


27% 
:, J—asin 
rese. À ê 


A(B)=b 


s 2H 
=) Rain 


> 
RO] 


et un calcul simple donne le résultat classique 


zasinG\ 


4 (6) — ab - 


z& Sin f 
À 


Considérons de même une base circulaire de 


rayon R (fig. 3) et cherchons sa directivité A (8). 


dans le plan XOZ. 


La fonction de répartition est égale à V2 Rx — a? 
pour o < x < 2R et nulle partout ailleurs. L'image 
de V2Rt —% est (©) 


o(p)— so eVARCpR) 


 T; étant la fonction de Bessel imaginaire du pre- 
mier ordre, on a donc 


97 p de 
I JR sin 5 27 : 
——_— ee À Ti =) > R sin $ 


2 « 
sin 


AKE)= mA 
Compte tenu de ce que Z() — jJ,(— j0, il vient 


UE (= sin s) 
A(B)=7R? ne ee 


2 


sin 8 


résultat bien connu. 
IL est à noter que la position de la base le long de 


l'axe des x n'intervient pas dans le résultat; une. 


translation de valeur + le long de l’axe des x se 
manifeste, en effet, dans l’image de la fonction de 
répartition par la présence d’un multiplicateur de 
la forme e-” dont le module est égal à l’unité avec 
- la valeur particulière que prend p dans les problèmes 
que nous étudions. 

Plus généralement, on peut remarquer, qu'étant 
_ donnée une base de fonction de répartition f (x) 
d'image 9 (p), toute fonction 


?(P)0(P) 

est l’image de la fonction de répartition F (x) d’une 
distribution ayant même directivité que la base 
précédemment envisagée si | 0 (p)| — 1 pour 


DR 
p=—) 7 sin 6. 


? 


() P. HUMBERT et McLacuLaN, Formulaire pour le calcul 
symbolique, p. 14 (Gauthier-Villars, Paris, 1941). 
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Soit, par exemple, un rectangle de hauteur b et 
de longueur a, comptée parallèlement à Ox:; l’image 
de sa fonction de répartition est 
Tr e //« 


P 


b 


Prenons Ü(p) — e / (x réel), qui est bien telle 
que |9 (p)| = r pour p = —j “sin. 
La fonction de répartition F(x) d’une distri- 


bution ayant même directivité que le rectangle dans 
le plan de tracé Ox a pour image 


œ (e 4 
à ne RER 
fie pa etes 2 
DEEE nee — h é-pa Ron 
P. ) 
Or, : 
LES 
er Er 
à E Jo (2 Var), 
donc 
bJ6(2 Vax) POUR 0 RQ, 
F(x) —= 


b ! Jo(2 Vax) — J(2 Vatz— a))} 


PORTER 


‘On obtient une distribution semi-infinie. 


2. Détermination de bases planes ayant des 
propriétés directives données. — Supposons que 


l’on se soit donné une base de fonction de répar- 
_tition f(x) et de directivité À (8) connues et que 
* l’on veuille déterminer la base plane ayant pour 


directivité À (8): (n entier positif), ce qui aura pour 
effet d’afliner le maximum principal (pris égal à 
l'unité) et de diminuer l’importance des maxima 


secondaires. 


Il suffira de chercher la fonction + (p) qui corres- 
pond à À (8) et qui n’est autre que la transformée de 
Laplace de f (x). D’après le théorème du produit de 
composition, la fonction de répartition de la base 
qui correspond à la directivité A (8): sera 


x 
J(&) 


x étant le symbole du produit de composition. 

Soit, par exemple, la base rectangulaire de lon- 
gueur a et de largeur b dont nous avons déterminé 
précédemment la directivité A (8) et cherchons la 
forme de la base possédant une directivité À (B)°. 
Sa fonction de répartition aura pour image 


b? fee ee } = 2? ( : 


D Fe 


LUE 


e—?2ap 9 e— «tp 
HR © — ——— |}. 
P? P? ) 


En remontant aux originaux, on voit que la 
fonction de répartition de la base cherchée est 


Bar pOur + 0% < 4, 
dÊRL IA, 


24 LT LH Pr 


4 
f(x) = | b(— x+92a) pour 


0 pour 
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La base a la forme d’un triangle isocèle de hau- 
teur b?a et dont la base, PAPReRE à Ox, a pour 
longueur 24. 

On aurait pu de même chercher la forme de la 
base ayant la directivité À (6. L'image de sa 
fonction de répartition est 


1e. ?4a n 
en (E) 
(1 5 


—= u( RU Le e- PAR - : e—2pa > 2 eue) 
P° D PE 


2 
ü 


P 
et l’on a 
b3 x? 
pour: o0<æ<4, 
2 
b3 
f(x) (22 +6ax—3 a) Pour a Tea, 
LT) = 2 
b3 
_ (æ?— 6ax + ga) pour 2a<x<34, 
0 pour 3QA<x<+o, 


La base est limitée par l’axe des x et un contour 
formé d’arcs de paraboles, ce contour étant symé- 


trique par rapport à la droite x — Sa. 


On aurait pu également se donner n bases de 
fonctions de répartition f; (x) et de directivités Ax (8) 
connues (= 12% .7n)etChercher là IOnC RON de 


répartition d’un groupement de directivité La. (B).. 


La fonction de répartition de la base cherchée est . 


[Le 


toujours par application du théorème du produit de, 


composition, 
Soient, par exemple, deux bases rectangulaires 
de largeur b et b' et de longueur a et a' (a > a). 
Les fonctions de répartition de ces bases_ ont 
respectivement pour images 
b nes , b' 
Le 1 


1—e P# 


et l’image de la fonction de répartition de la base 
cherchée s'écrit 


[x — e—P4_ pe—pa' LE e—pla+a as 
P° 


bb" 


Cette fonction f (x) s’écrit donc 


bb'æ pour DIT TION 


() bb'a pour HELD EROR 
& he Ê ==} 
: bb'(— x + a+ a) pour a'Lzx<a+tt, 


0 pour a+ dd <LxL+ 


et la base cherchée a la forme d’un trapèze de tracé 
facile. 


3. Généralisation. — Plus Dents fenont consi- 
dérons un groupement constitué par n bases na 
rectangulaires identiques À, équidistantes et dispo- 
sées comme l'indique la figure ds tous leurs points 


vibrants en phase. 4 


Fig. 4. 4 


Les divers rectangles. se déduisant chacun du à 


précédent par une translation parallèle à OX de 


valeur 1, l’image de la fonction de répartition de ce 
groupement relative au plan de trace OX est 


rene Cane 
e(p)=b} = To Me 
re TT se | 
PA 
1— e—Pa | 
= b————[r1+e-rl+,.,.+ ep]. 


La directivité À (8) du groupement dans le plan 2 
de trace OX en résulte immédiatement; ïl vient 


AGE) = |, 
soit 
ee (ESS) Le (= rn sin PE 
ANS MELUÉ 
À ee 


Toujours dans le but d'obtenir une meilleure : 
directivité, supposons que l’on veuille déterminer la 
structure du groupement qui aura pour directivité 
une puissance donnée, m, de celle du groupement 
initial dans le plan de trace OX. L'image de sa 
fonction de répartition sera 


d(p) = p(p}— om)" 


xX[i1+e-Pi+,,.+e-pti-ti]n, 


Le premier terme du produit qui figure dans 4 
l'expression de © (p) n’est autre que l’image de la 


fonction de répartition d’une base A’ ayant, dans le 
plan de trace OX, une directivité égale à la puis- 
sance m de celle de l’un des rectangles qui cons- 
tituent le groupement initial. 

Le second terme est de la forme 


k=m{(n—1) 


> ax er PA, 
k=—0 


les «x étant les coefficients du développement du 
polynome 


(LH 224. ,+ mit), 


On voit ainsi que la nouvelle base sera constituée 
par 1 + m(n—1:) éléments identiques à A’ qui se 
déduisent chacun du précédent, comme dans la base 
initiale, par une translation de valeur ! parallèle 

à OX, leurs poids devant être pris Fu aux coeffi- 
{cents correspondants «x. 

* A titre d'exemple, considérons une Da à deux 
Lente À et cherchons la structure du groupement 
qui aura dans le plan de trace OX, une directivité 
égale au carré de celle du groupement initial. 

La fonction de use du groupement cherché 
aura pour image 


Le e= pas 
b? (==) ({i+2e-Pl+ e-2pl) 


et il apparaît que ce groupement sera constitué par 
trois triangles isocèles identiques A’, de hauteur b?a 
_et de base 24, chacun d’eux se déduisant du précédent 
par une translation de valeur / parallèle à OX, les 
poids respectifs de ces bases partielles étant 1,2 et r. 
On peut, pour tenir compte de ces poids, imaginer 
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une base à quatre triangles ayant la 


structure 
figurée en 5. 


Fig.:.5. 


De ce résultat, il ressort, qu'étant donnée une 
base à deux éléments, identique à celle que nous 
venons de considérer, on peut, en prenant une 
valeur de m suffisamment élevée, déterminer un 
groupement linéaire comportant un plus grand 
nombre d’éléments et ayant une directivité aussi 
satisfaisante que l’on veut. 


Ce dernier résultat est à rapprocher de celui 
obtenu par Shelkunoff dans son article précité. 


Manuscrit reçu le 26 mars 1948. 
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VARIATION THERMIQUE DU CHAMP COERCITIF DU COBALT DIVISÉ 


Par L. WEIL, S. MARFOURE et F. BERTAUT. 
Laboratoire d’Électrostatique et de Physique du Métal. 


Sommaire. — On a étudié entre la température de l’air liquide et 3000 la variation du champ 
coercitif de cobalt divisé, obtenu par réduction d’oxalate, de formiate et d’hydroxyde et comprimé 
en forme de bâtonnets. On a déterminé par un examen aux rayons X la structure et, en particulier, 
la proportion de phase hexagonale et cubique de chaque échantillon. 

Pour tous les échantillons à forte proportion de phase hexagonale, le champ coercitif passe par un 
minimum au voisinage de 2000 et a alors une valeur quatre fois plus petite que dans l'air liquide. 
Ce résultat est conforme à la théorie de M. Néel qui attribue le champ coercitif des poudres en partie 
à l’anisotropie magnétique qui, pour le cobalt hexagonal, s’annule au voisinage de 200° et croît quand 


la température décroît. 


Nous avons montré, dans une publication anté- 
rieure [r] que le champ coercitif du nickel divisé, 
obtenu par réduction chimique et aggloméré vérifiait, 
d’une manière tout au moins qualitative, la théorie 
du champ coercitif des ferromagnétiques à grain fin 
développée par M. Néel [2, 3]. 

D'après cette théorie, le champ coercitif résulte, 
. d’une part de l’anisotropie de forme des grains, d'autre 
part de l’anisotropie magnétique de la substance. Les 


valeurs qu’on peut s'attendre à trouver pour le 
cobalt hexagonal en grains fins isolés sont respecti- 
vement de l’ordre de H} — 0,26. T.ÿ — 1170 gauss 


K ; à 
et H}— 0,95 + — 2900 gauss (1) à température 


ordinaire, en prenant pour K la valeur 4,1.10% (4,6). 
Pour le cobalt cubique, H{ serait le même, tandis 


() Formule inédite communiquée par M. Néel. 
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ES 6 HER 
ue H° serait donné par la formule 0,64 où K est 
q P # 


cette fois la constante d’anisotropie de la phase 
cubique qu’on peut situer, par extrapolation des 
valeurs obtenues pour les alliages Fe-Co et Ni-Co 
cubiques, au voisinage de 1.106, ce qui donne 
H} = 500 gauss. Nous avons montré, dans le cas 
du fer et d’un ferro-cobalt à 30 pour 100 [5], dont 
l’anisotropie magnétique est faible, que le champ 
coercitif des comprimés est toujours bien inférieur 
à la valeur de H!/ calculée. Il y a tout lieu de penser 
que le champ coercitif dû à l’anisotropie magnétique 
est également réduit lorsque l’on rapproche les 
grains au Cours de l’agglomération. 

Mais comment distinguer l’effet de la forme des 
grains de l’effet de l’anisotropie magnétique ? H/ varie 
lentement avec la température et s’annule au point 
de Curie, soit à 11000 pour le cobalt. H”* doit, au 
contraire, varier considérablement : on sait, en 
effet [6], que X augmente d’environ 5o pour ro0 
entre 12° et — 1839 et s’annule pour changer de 
signe aux environs de 2600; à 390°, K reprend, 
par exemple au signe près, la même valeur qu’à 1500. 
En raison de la superposition de H/ et H”, on doit 
donc s’attendre à trouver un minimum non nul du 
champ coercitif au voisinage de 2600. 


Les échantillons. — Le principe de notre méthode 
de préparation dû à MM. Néel, Aubry et à l’un de 
nous est connu [8, 1] : on décompose et réduit un 
oxalate, un formiate ou un hydroxyde cobalteux 
dans l'hydrogène à basse température. La poudre 
obtenue est comprimée en bâtonnets. Ceux-ci sont 
recuits à une température un peu supérieure à celle 
qui est atteinte au cours de l’étude de la variation 
thermique. 

Nous avons ps élément sur des échantillons 
de cobalt avec 5 pour 100 de zinc. L’oxalate utilisé 
pour les obtenir a été préparé par Préc RUE d'une 
solution de sel de cobalt additionné de 5 pour 100 
de sel de zinc. Nous recherchions l'influence du zinc 
sur la structure. 


La structure du cobalt. — On connaît le cobalt 
sous deux formes cristallines, l’une hexagonale (+), 
stable au-dessous de /4oo°, l’autre cubique (y), 
stable à haute température. La transformation & > y 
s'accompagne d’un faible effet thermique (62 cal par 
atome-gramme [13]) et est lente; il est relativement 
facile d'obtenir du cobalt cubique à la température 
ordinaire par un refroidissement rapide par exemple. 
De plus, la plupart des impuretés abaissent le point 
de transformation. On a montré [o| qu'il suffisait 
de 0,08 pour 100 de carbone pour l’abaisser à 14920. 
Edwards et Lipson [ro] ont même mis en évidence 
que, pour du cobalt pulvérisé mécaniquement en 
grains de moins de 5o p, il était impossible, même 
par des recuits prolongés — au-dessous de la tempé- 
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rature de recristallisation toutefois — de dépasser « 
une teneur de 5o pour 100 en phase +. Enfin, d’autres 4 
auteurs ont signalé que, d’une façon générale, les M 
tensions internes dues à l’écrouissage [11] ou les M 
inclusions [12] stabilisaient la phase cubique et " 
pouvaient empêcher un rétablissement, même partiel, : 
de la phase hexagonale au cours du refroidissement. « 

II était donc indispensable de s’assurer de la struc- M 
ture des échantillons étudiés pour pouvoir inter- » 
préter leurs propriétés magnétiques. ; 


La structure des échantillons. — Nous avons - 
fait des diagrammes de Debye-Scherrer des divers 
échantillons étudiés. En employant la méthode de 
focalisation de J. C. Brentano, nous avons déduit la … 
proportion de cobalt hexagonal de la mesure du 
rapport des intensités de deux raies; l’une (ro11) 
est spécifique de la forme hexagonale; l’autre corres- 


pond soit à (0002) dans la forme hexagonale, soit | 


à (x11) dans la forme cubique. 

On sait [14] que l’on passe de la phase hexagonale « 
compacte à la phase cubique par un simple glis- : 
sement, parallèle à un axe binaire, 


gueur NE des plans (ooo1) qui deviennent des 


plans (111) du cube. Or Edwards, Lipson et 
Wilson [15] ont montré que, dans la phase hexa- 
gonale, des séquences de plans (ooo1) peuvent être 
décalées l’une par raprort à l’autre par un glis- 
sement semblable. Tout se passe donc pour certains 
plans réticulaires tels que (ooo1) et, d’une façon 
générale, tous les plans communs au cube et à la 
phase hexagonale, comme si l’on avait affaire à un 
grain d’une certaine dimension, alors que, pour les 


. autres, les largeurs des raies de diffraction corres-. « 


pondent à des feuillets d'épaisseur bien plus petite. 

On peut ainsi, par l'étude des largeurs de diverses 
raies, évaluer la largeur et l’épaisseur moyenne des 
feuillets. Nous avons trouvé, pour les échantiilons 1 
à 3 contenant essentiellement la phase hexagonale, 
des épaisseurs de 30 À pour des largeurs de 250 À 
environ, la hauteur de l’empilement étant de l’ordre 
de 250 À. 

Le tableau: ci-dessous résume les observations 
faites sur les divers échantillons et leur mode de 
préparation. 


Variation thermique du champ coerciti. — 
Le Tableau IT ci-contre donne les valeurs du champ 
coercitif pour les dix échantillons à diverses tempé- 
ratures dans l’ordre où elles ont été obtenues par 
la méthode antérieurement décrite [1]. Nous avons 
reporté les résultats sur les figures 1 et 2. 

Pour les échantillons 1 à 3 et 7 à 10, d’une façon 
très marquée, pour les échantillons 4 et 5 d’une 
façon moins prononcée, le champ coercitif passe 
par un minimum au voisinage de 2000. Ce minimum 
n'existe pas pour le cobalt cubique. 


d’une lon- = 


"Rat PNee à 
mi 


à 


_ VARIATION THERMIQUE DU CHAMP COERCITIF DU COBALT DIVISÉ 203 


TABLEAU [. 
Échantillon. Préparation, Recuit. Structure, 
Es Réduction d’oxalate pur, Sous vide, 75 min. à 3000 


1/2 hexagonal 
50 min à 5400 


DES » D 48 heures  »: 2/3 » 
ASE ENS Oe » e » 15 jours » 2/3 » 
AR .. Réduction d’oxalate de Co » 2 h 45 à 5000 Peu d’hexagonal 
+ 5 pour 100 Zn à 3500 
RO Id., à 3400 » 2 h 15 à 360° RAS 
Ge e Id., à 3000 » à 300° Cubique ; pas d’hexagonal 
à décelable 
rare Réduction de formiate pur, Dans l’hydrogène, Peu de cubique 
5o min à 3200 1 heure à 3500 
REREARE _ Réduction d’hydroxyde cobalteux, Id., 2 heures à 3500 » (Ai 
5o min à 3200 
D [d.;:à 3500 °° » À » 
AO ESS Id., à 4502 | » » 
._ (*) Avant recuit cet échantillon était hexagonal pur. 
TABLEAU Il. 
af: li Pie 3. 4. 5. 
CR NE ie << 2 RS CR A CR SE 
t H, t H, as H, £ AH. t A, 
—190 640 —190 706 —190 663 —190 178 T0 100 
— 79 433 — 79 MAT 13 340 102 104 — 79 296 
19 322 13 322 75 265 D. 100 15 245 
66 265 71 235 F0: 005 58 98 59 235 
104 246 [19 203 TADSS0 7D 98 71 234,4 
128 230 149 182 162. 250 88 98 90 234 
184 ptet 189 172 CDS MATOS 100 99 122 234; D 
232 220 238 179 LORS AIO ED 118 IOI 131 235 
284 235 208 200 202 183 140 102,5 157 236 
208 235 48 181 252 186 128 104 228 237 
; 218 172 230 185 208 109 102-7209 4 
164 179 202 183 258 TES SI 234 
170 191 310 119,9 10 244 
123 220 350 117 
6 7 8. 9 10 
Re ns DRRER EAN PRE ES 5 EUR a rene ATEN ES 
3 Û IT, . HE t H, t H t I, 
—190 790 319 197,9 192 250 13 239 13 208 
15 6) 270 190 Ho 197 ; 5 63 192 65 0) 
60. 554 211 183 108. Fb3,0 126 128,2 123 124 
82 548 171 186,9 TPE 125 166 11270 186 109,7 
102 548 121 206 249 130,9 214 109,7 220 106,5 
149 544 73 238 309 163.,5 270 120,0 265 116,3 
210 535 13 301,9 1 247 364 149 314 128,5 
290 518 =-190 565 —190 450,5 13 242 13 204 
175 540 mL) 408 #10) Dr run 510 190 12 
105 548 13 300 A T) 397,9 — 79... 298 
HO 007 
Conclusion. — Le minimum de champ coercilif théorie. Il met bien en évidence la double origine du 


\ 


observé pour tous les échantillons à forte propor- champ coercitif des poudres ferro-magnétiques 1e 
tion de phase hexagonale est en accord avec la minimum se produit pres du point où la constante 
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d’anisotropie s’annule; ce minimum n’est pas nul 
parce que l’anisotropie de forme contribue au champ 
coercitif. 


800} Ï + 
He 6 
700 e 
600 
500 LR 
400 
DS 
200 LR | 
; ee 
Î 
100 ; PRES ESS SR Canard EEE 
- 100 100 200 S00E TR 
Fig. 1. 


Pour l'échantillon 6, purement cubique aux 
rayons X, le minimum est remplacé par un point 
d'inflexion, dû peut-être à des traces de phase y 


. non décelées; d’ailleurs, pour le cobalt cubique, 


on ne sait rien de la variation thermique de K. 


Le minimum a toujours lieu à une température 
nettement inférieure à 2600, alors que la superpo- 
sition de H/, décroissant régulièrement avec la 
température, à H° devrait plutôt le décaler vers 
les températures élevées. Nous pensons que la 
variation thermique de K pour nos échantillons 
est différente de celle observée pour le cobalt massif. 
Tout se passe comme si la température où K s’annule 
était abaissée par les impuretés (carbone, oxygène, 
zinc, etc.), la structure feuilletée ou même simplement 
la finesse du grain. Ce résultat n’a rien de surpre- 
nant : le point de transformation y -> € est lui-même 
abaissé par ces causes ainsi que le rappelle la persis- 
tance tout au moins partielle, à température ordi- 
naire, de cobalt cubique dans nos échantillons, 
cépendant que le minimum de H, se situe d'autant 


plus bas que les re Sont plus os descen- . 
dant, par exemple, à 1000 pour les éprouvettes 4. 
et 5, contenant beaucoup d’impuretés et de y. 

Un véritable décalage dans la variation ther- | 
mique de K permet également de comprendre 4 
l’augmentation extrêmement importante de H. entre , 
la température ambiante et celle de l’air liquide : \ 
à une augmentation de 50 pour 100 de X dans cet | 
intervalle comme on l’a observée pour le métal massif 
ne saurait correspondre un doublement ou un tri- 0 
plement de champ coercitif comme nous l'avons à 
observé pour nos échantillons. 


Hele 
gausS 


500 


200 


100 


-/00 0) 100 200 300 t= 


. Reste à savoir pourquoi les valeurs de H. sont « 
plus faibles que celles qu'on s'attendait à trouver . 
(moins de 4oo au lieu de 2900 à température 
ambiante !) et pourquoi les valeurs pour la phase y - 
sont du même ordre, ou même plus grandes que 
celles pour la phase hexagonale. Les valeurs obser- 
vées sont inférieures aux valeurs calculées en raison 
du rapprochement des grains comme nous l'avons 
déjà signalé plus haut. Quant au faible champ « 
coercitif du cobalt hexagonal, il trouve sans doute « 
son explication dans une micrographie électronique 

de J. T. Me Cartney et R. B. Anderson : ces auteurs # 
ont montré que les grains de cobalt hexagonal étaient 
de larges lames parallèles au plan (ooo1), groupés 
parallèlement et non orientés au hasard comme on 
aurait pu le penser a priori. Le fort champ démagné- “ 
tisant parallèle à l’axe sénaire empêche, tout au 
moins en partie, celui-ci de jouer son rôle de direc- * 
tion de facile aimantation. On perd donc presque M 
complètement le bénéfice de la forte anisotropie 
magnétique qui lui était associée, puisque, dans le . 
plan de base, le cobalt est sensiblement isotrope [4]. 


Manuscrit reçu le 12 avril 1948. 


2 


ces < 


= ÉTAT DES CONTRAINTES AUTOUR DES ACCROCHAGES 


FRS ; BIBLIOGRAPHIE. 


| SETIL: do MARFOURE, J. de Physique, 1947, 8° série, 
ÿ P: 990. 

[2] L. Née, C. R. Acad. Sc., 1947, 224, p. 1488. 

- [3] L. Néez, C. R. Acad. Sc. 1947, 224, p. 1550. 

[41 R. Becker et W. DôRING, Ferromagnetismus, Springer, 
le Berlin, 1930. 

[5] L. Werz, C. R. Acad. Sc., 1947, 225, P. 220. 

[6] K.: Honpa et H. Masumoro, Sci. Rep. Tokohu Uni. 

3 1991,20,; p: 323. 
[7] R. FORRER, R. BAFRIE et P. FOURNIER, J. de Physique, 
| | 1945, 8e série, 6, p. 50. 
_ [8] Brevet français, Chambéry, 9 avril 1942. 


ne Aufbau der Zweistofflegierungen, Springer, 

1936. 

[ro] O. S. Enwarps et H. Lipson, Tech. J. of Inst. of 
Metals, 1943, p. 177. 

[11] G. WASsERMANN, Metallwirtschaft, 1932, 41, p. 6r. 

[121 R. Fricxe et H. MüLcer, Naturwissenschaften, 1942, 
30, p. 439. 

[131 U. DEHLINGER, Z. f. Physik, 1931, 68, p. 535. 

[14] BraGG, Crystalline State, p. 144 (Bell, London, 1933). 
[15] O. S. Epwanps, H. LrpsoN et A. J. C. Wizson, Nature, 
1941, 148, p. 165. 

[16] J. T. Mc CarrNey et R. B. ANDERSON, J. appl. Physics, 

1947, 18, p. 902. 


907 


_ [91 HosxrmaTo, Kinzoky no Kenkyu, 1932, 9, p. 57; d’après 


| ë k : 
LE ED TENTE EN ERENE ONE LEE CE DENT A TES PEN RESTE EP RER EE RG  P 


LE JOURNAL DE PHYSIQUE ET LE RADIUM. 


SÉRIE VIII, TOME IX, JUIN 1948. 


SUR L'ÉTAT DES CONTRAINTES AUTOUR DES ACCROCHAGES 
DE LA THÉORIE DE LA PLASTICITÉ DES MÉTAUX 


Par PIERRE LAURENT, 
Chargé de Recherches du C. N. R.S. 


Sommaire. — Rappel des propriétés plastiques des cristaux et des théories de Smekal, Becker et Taylor. 
En admettant que la déformation plastique est due à la formation, à la propagation et à la disparition 
de défauts du réseau cristallin ou accrochages, on obtient un bon accord avec les résultats expéri- 
mentaux. Cette théorie a été précédemment développée [3] en se donnant, «a priori, la répartition des 
contraintes internes autour d’un accrochage; dans cette note, il est montré que l'hypothèse faite est en 


accord avec les conditions de l'équilibre élastique. 


. Les alliages étant constitués de cristaux, leurs 
propriétés dérivent des propriétés de ceux-ci et de 
leurs interactions. L'étude de la plasticité des 
cristaux uniques est ainsi un premier stade pour la 
compréhension des caractéristiques mécaniques. 


_ Plasticité des cristaux uniques {[1]. — Tant que 
les forces appliquées à un cristal sont faibles, les 
déformations suivent la loi de Hooke et le cristal 
reprend sa forme initiale quand on supprime les 
forces extérieures. 
_ L'expérience montre que si l'importance des 
forces appliquées augmente, le cristal ne reprend 
plus sa forme initiale; on dit qu’il a subi des défor- 
mations plastiques. Celles-ci peuvent être de trois 
types : 

a. déformation par translation; 

b. maclage mécanique; 

c. kinking ou pliage en genou. 


Le premier mode de déformation est extrêmement 
général et de beaucoup le plus important; le maclage 
ne s’observe que sur certains types de cristaux, 
quand au kinking il semble encore exceptionnel. 


Les déformations par translation sont dues à des 
glissements d’une partie du cristal par rapport à 
l’autre le long d’un plan T, cristallographiquement 
défini (plan de translation); dans ce plan, la direc- 
tion de glissement est parallèle à une direction: { 
également cristallographiquement déterminée (direc- 
tion de translation). Si d est la distance de deux 
atomes voisins dans la direction f, la grandeur de la 
translation est toujours un multiple de d. 

Si l’on enregistre des diagrammes de traction 
pour différentes orientations de l’axe de traction 
par rapport au réseau du cristal, on obtient des 
courbes très diverses et des limites élastiques très 
différentes. Mais on peut constater que le cisail- 
lement agissant sur le plan de translation actif au 
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moment dé l'apparition du premier glissement plas- 
tique est indépendant de l'orientation du réseau; 

telle est la loi du cisaillement critique de Schmid. 
Si, au lieu d'envisager comme on le fait habituel- 
lement la variation de l'effort de traction en fonc- 
tion de l'allongement, on considère la variation du 
cisaillement agissant dans la direction de translation 
en fonction de la valeur du glissement, on obtient 
ainsi la «courbe de consolidation » du métal indépen- 
dante de l'orientation du cristal; elle représente 
une caractéristique du métal considéré. 


Théories de la plasticité cristalline. — Pour 
expliquer la résistance des cristaux, on peut chercher 
à calculer le travail nécessaire pour rompre les 
liaisons interatomiques; ce travail a été entrepris 
par Born. Toutefois, les valeurs obtenues étant 
mille fois plus élevées que les résultats expérimentaux 
et la validité du raisonnement de Born ne pouvant 
être mise en doute, il est nécessaire d'expliquer les 
faibles résistances observées en pratique. Les 
recherches entreprises dans cette voie ont donné 
naissance à diverses théories parmi lesquelles nous 
citerons : 


1° La théorie de Smekal, qui admet l’existence de 
microfissures dans le cristal. Ces fissures provoquent 
des concentrations d'efforts entraînant localement 
le dépassement de la résistance théorique. Toutefois, 
étant donné les dimensions admissibles pour ces 
microfissures, on ne peut combler le fossé entre la 
résistance pratique et la résistance théorique de 
Born. | 


20 La théorie de Taylor, qui admet l'existence dans 
le réseau d’atomes dont les liaisons avec leurs 
voisins sont anormales; ceci constitue ce que l’on 
appelle une dislocation. Ces dislocations provien- 
draient d’un vice congénital du cristal. Sous l’effet 
d’une force extérieure relativement faible, cette 
dislocation peut se déplacer, ce qui entraîne un 
glissement relatif d’une maille d’une partie du cristal 
par rapport à l’autre. Cette théorie explique bien 
le mécanisme du glissement et l’existence du cisail- 
lement critique, mais ne permet pas de comprendre 
l’influence de la température, ni de la vitesse de 
déformation. 


39 La théorie de Becker qui fait appel non pas à 
une hétérogénéité de structure comme celle de 
Taylor, mais à l’hétérogénéité des énergies vibra- 
toires des atomes. Si un atome possède une grande 
énergie, l'amplitude de ces vibrations est également 
grande par rapport à la distance réticulaire; aussi, 
sufiira-t-il d'appliquer une faible force extérieure 
pour que cet atome quitte sa position d'équilibre 
normal et prenne la place d’un atome voisin. Si ce 
fait se répète tout le long d’un plan de glissement, 
il en résulte, pour l’ensemble du cristal, un glis- 
sement d’une maille, 
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Une synthèse de ces théories a été entreprise 4 
successivement par Orowan [4],par W.G. Burgers [2] « 
et par Kochendorfer [5]; elle a permis d’expliquer « 
l'effet de la température sur le cisaillement critique . 


E, 


et sur la pente à l’origine des courbes de consolidation 


aux températures moyennes, ainsi que de calculer la 


courbe de consolidation de l'aluminium. Ces théories 


sont restées impuissantes pour retrouver les courbes 


de consolidation des cristaux hexagonaux et pour = 


expliquer le fluage, la plasticité au zéro absolu, 


la constance du cisaillement critique et du coeffi- 
cient de consolidation au voisinage du point de … 


fusion, ainsi que la restauration et la stabilité de « 


l’état déformé. La raison de ces difficultés est due, 


à notre avis, à ce que l’on a adopté comme un 
postulat des propriétés les dislocations proposées 
par Taylor. 


Théorie des accrochages. — Dans cette théorie 
on admet, comme dans-la théorie de Taylor que le « 
glissement plastique résulte de la propagation à . 
travers le cristal des défauts du réseau que nous = 
appellerons des accrochages pour les distinguer des " 


dislocations de Taylor; nous nous contenterons de 


résumer ici les bases de cette théorie [3]. Ces accro- 
chages se forment à la limite du bloc mosaïque sous 


l’action combinée de l'agitation thermique et du 
cisaillement agissant sur le cristal; ils se propagent 
dans le plan de glissement parallèlement à la direc- 
tion de glissement et viennent s’accumuler à l’autre 
extrémité du bloc mosaïque. Ces accrochages créent 
autour d'eux une perturbation du réseau, c’est-à-dire 
un champ de tensions internes qui, en s’opposant 
au cisaillement extérieur, est à l’origine de la conso- 
lidation du cristal. Au fur et à mesure que la défor- 
mation augmente, le nombre des accrochages accu- 
mulés croît, mais un certain nombre de ceux-ci est 


détruit par l'agitation thermique (restauration). La 
disparition d’un accrochage entraîne la disparition 


du champ des contraintes qu’il créait, mais ne 
modifie pas le glissement plastique que la propa- 
gation de cet accrochage a provoqué. 


Contraintes autour d'un accrochage. —_ En 
développant cette théorie, on admet que les ten- 


sions internes créées par un accrochage s’annulent « 


pratiquement à une certaine distance du centre de 


l'accrochage et que le cisaillement dans le plan de M 


glissement s'oppose au cisaillement extérieur sur les 
atomes situés du côté d'où cet accrochage vient et 
s'ajoute au cisaillement extérieur sur les atomes 
situés du côté où il se dirige. 

Plaçons-nous pour simplifier dans le cas -d’un 
cristal à deux dimensions et prenons comme origine 
des angles la démi-droite que vient de parcourir 
l'accrochage. Les hypothèses précédentes peuvent 
être satisfaites en admettant qu’en un point quel- 
conque le cisaillement + créé par l’accrochage est 


K cos 


LAS 


Gr 


1 


où r et 0 sont les coordonnées polaires du point 
_ considéré en prenant comme origine la ligne de 
glissement, K et «x des constantes. Nous suppo- 
serons, de plus, l’accrochage non pas dans un bloc 
mosaïque, mais dans un plan infini, c’est-à-dire que 
_ les composantes du Don de contraintes internes 0, 


_ so et Tr s’annulent avec - ; 


Que peut-on dire du champ des déformations 
autour de l’accrochage ? J. M. Burgers [6] fait 
remarquer que, pour une dislocation, si à la dis- 
_ tance r la composante du déplacement parallèle 
à un rayon polaire est u pour 9 —o, il a la 
valeur u + d pour 0 — 2 x (d étant le paramètre du 
réseau dans la direction de glissement); la compo- 
sante u du déplacement est ainsi une fonction 
multiforme de l’angle polaire. Au contraire, la compo- 
sante v du déplacement normal au rayon polaire 
est une fonction uniforme. Burgers étudie alors un 
état de contrainte élastique satisfaisant à ces condi- 
tions. Ce point de vue est incompatible avec les 
propriétés attribuées aux accrochages. En effet, 
_ lorsque dans un état de tensions élastiques on fait 
disparaître les contraintes en tous les points, les 
déformations élastiques disparaissent ipso facto. 
Or, nous avons, au contraire, supposé que, lorsqu'un 
_accrochage disparaît (restauration) et que le champ 
de tensions s’annule, le glissement de la partie 
supérieure du cristal par rapport à la partie infé- 


 rieure reste égal à d. Donc le déplacément u, fonc- 


tion multiforme de 0 ne peut alors être considéré 
- comme un déplacement élastique. Donc, pour 0 — 27, 


le déplacement total se compose d’un déplacement 


- non élastique d et d’une partie élastique (que nous 
désignerons par u fonction uniforme de 0; le 
. déplacement non élastique est, au contraire, une 
fonction multiforme de 0). 

Le déplacement d résulte du mécanisme propre 
- de propagation de l’accrochage, nous ne le considé- 
rerons pas ici. Le champ de contraintes élastiques 
autour d’un accrochage doit satisfaire aux condi- 
_tions suivantes : 


10 le cisaillement le long du rayon polaire est 


K cos®. 
ra? 


T 


20 les CORPORANES des contraintes s’annulent 


avec + = 


SA _. déplacements élastiques sont des fonctions 


à . I 
uniformes de 0 et s’annulent avec = 


En coordonnées polaires, les équations de l’équi- : 


libre élastique sont, en l’absence de forces mas- 
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siques [7] : 
d O7 
ra 0; (2) 
Ole) 
0 + e 
PT go Te T'es 0 (3) 


Par suite de Laine (1) faite sur 7, l’équa- 
tion (3) s'écrit 
V6 __ Kcosû 
Ha de 


(a —2), 


ce qui donne 
K sin > 


ds 
Sp— (a12) ap 0 0 (4) 
où f(r) est une fonction de r seul encore indéter- 
minée. Portons cette valeur de 55 dans l’équation (2); 
il vient après intégration 
116) sd k sin 
7 Î 


r% 


y , 
où g (0) est une fonction de 0 seul encore indéter- 
minée. L'état des contraintes est ainsi 

in 0 
ne ee Pet) 


RAA COS 
ne 


Ù 


K sin 3 
tr 

Pour déterminer f (r) et g (0), il faudrait écrire que: 
les fonctions 5,, on et 7 satisfont à l’équation de 
comptabilité; il est plus simple de calculer direc- 
tement les déplacements u et v. 


Détermination des déplacements autour d'un 
accrochage. — On sait que les déformations &,, 
ep et y sont liées aux contraintes par les relations 


Fe; = .5r— V5 Eh = 54 — vor, Ey=2(1+ v}, 
où E et » sont le module de Young et le coefficient 
de Poisson, tandis que les déformations sont reliées 
aux déplacements u et v par les relations 
1 du Fe 
r 08 ? 

1 Ou 06.0 | 


Aves du sr 
nant 
(6) 


té n08 or re 
Portons les valeurs (5) dans ces équations; on trouve 
alors 
u=] A0 dr + g(8)Lnr — vf(r) 


” . ( Fo D y y A 
MORE Are tar NS (7) 
ra—1 CEA | 


où À (0) est une fonction indéterminée encore de 0 
seul. Remarquons que u ne peut tendre vers zéro 


avec = pour toutes les valeurs de 0 que si 
h() = &(8) EG) 


210 


et si et f (r) tendent au zéro avec =. La deuxième 


équation (6) s’écrit, dans ces conditions, 


+0 fr 0 far) + ç, (8) 


où L(r) est une fonction encore indéterminée de r 
“ I : AE 
seul, s’annulant avec -- Mais, puisque v doit être 


une fonction uniforme de 0, il faut que l’on ait 


A re AGAIN SE 
2 es rt 


cette équation différentielle du type de Fuchs a 


pour solution une fonction de la forme © + br 


(où a et b sont des constantes). D’après les condi- 
tions imposées à f(r), on a b — o.. Portons les 
valeurs (7) et (8) dans la troisième équation (6) 
en tenant compte de (9); on obtient une expression 
| : K cos 0 
18 


contenant des termes en et des termes 


dépendant de r seulement; elle ne peut être satis- 
faite que si les deux termes sont séparément nuls, 
c’est-à-dire à la fois 


K cos 0 


rx 


d 2 Le 
PUR) T0; (a?— 2-3) =0; 


(10) 


la première équation n’admet que la solution ! — 0, 
PS I Fa . 
qui s’annule avec = : la deuxième donne soit x = 3, 


soit &« = — 1. La solution « — — 1 ne satisfait pas 
au problème. On a finalement, pour le champ de 
contraintes : 


a K sin a K sin 0 
TS ro Gps ce INR 
7? 7 r2 Fr u) 
_ Kcosû ( 
es 


Remarquons que, dans le raisonnement précédent, 
nous avons admis que l’accrochage était situé dans 
un plan indéfini. Notre calcul n’est donc plus 
valable en toute rigueur pour les accrochages situés 
_dans les blocs mosaïques des bords du cristal (. 
D'autre part, nous avons implicitement assimilé le 
cristal au milieu continu en écrivant les équa- 
tions (2), (3) et (6); ce fait est certainement inexact 
dans le domaine voisin du centre de l’accrochage. 

Dans la suite du raisonnement, nous suppose- 
rons a — o, cette hypothèse ne se justifie d’ailleurs 
que par la simplification des calculs qu’elle entraîne. 


(:) Ge cas pourrait se traiter par un calcul analogue, au lieu 
d'admettre que le champ de tension s’accumulant avec, 
il suffirait de considérer qu'il disparaît sur la frontière du 
cristal, a" 
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(9): 


Calcul de la consolidation atomique. — Les . 
résultats précédents nous conduisent à reprendre le « 
calcul de la consolidation atomique fait précé- # 
demment [3]. Ainsi que nous l'avons rappelé, les … 


accrochages accumulés créent, au point de forma- 


tion des accrochages, un cisaillement 7; qui s'oppose 
au cisaillement extérieur 7. Si 7, est le cisaillement 
nécessaire à la formation des accrochages, on aura 
T—T; = TO. à (192) 

r; est, par définition, la consolidation résultant de … 
la déformation déjà imposée au cristal. 4 
Ce cisaillement antagoniste qui s'exerce au point M 
de formation À des accrochages (fig. 1), peut se « 


Limite des blocs-mosaique Limite des blocs-mosaïque 


Direction hu glissem* 


— | — <— 
+ Accrochage direct =» Courant des accrochages directs 
_ F INVEISE ——. > inverses 

Fig. 1. — Accumulation des accrochages 


aux limites des blocs mosaïques. 


décomposer en deux parties : 


10 Cisaillement dû 
en a et b;; 


aux accrochages accumulés 


20 Cisaillement dû aux accrochages accumulés 


en B. 


Avant d'évaluer séparément ces deux termes, 
rappelons qu'il existe deux sortes d’accrochages, se 
déplaçant en sens inverse sous l’effet du cisaillement 
extérieur, les accrochages directs et les accrochages 
inverses. 


19 Cisaillement antagoniste dû aux accrochages 
accumulés en & et b;j. — Soient N4 le nombre des 
accrochages directs et N; celui des accrochages 
inverses accumulés en a; et b;. D’après ce qui précède 


les accrochages directs accumulés en @ créent, # 


en À, un cisaillement positif et ceux accumulés « 
en b;, un cisaillement négatif; inversement, les « 
accrochages inverses accumulés en à créent, en À, 
un cisaillement négatif et ceux accumulés en b; un 
cisaillement positif. Par suite, le cisaillement anta- 
goniste résultant de ces divers accrochages est 


m D ! 
RUE I I 
MEN AT RDS | 
“ 1 


m et p étant les nombres de blocs mosaïques situés 
à droite et à gauche de la frontière considérée, 


(13) 


10 


N° 6.. | | ÉTAT DES CONTRAINTES AUTOUR DES ACCROCHAGES 911 


| 


IL la dimension du bloc mosaïque. Les deux sommes 
‘entre crochets étant rapidement convergentes, on a 
paiement 


C}= (Na N) PE Vers 


(14) 


Ch est ainsi indépendant des be m et p, 
c’est-à-dire des dimensions du cristal. Pour les blocs 
| mosaïques situés au voisinage des bords du cristal, 
le calcul précédent n’est pas valable, mais nous 
n’envisagerons pas ici l'influence des dimensions du 
cristal. 


2° Cisaillement antagoniste dû aux accrochages 
accumulés au voisinage de B. — Pour calculer (r;);, 
 déterminons d’abord le cisaillement créé en A 
(fig. 2) dans la direction de glissement par un accro- 


A .Pont de formation des accrochages 


B. , d'accumulal! , no 


OI0cS- moSaque 


_—_— Courant des accrochages directs 


RS m »  lnverses 


mines = = = im 


Fig. 2. — Disposition des accrochages 
à la limite de deux blocs mosaïques. 


chage d’abscisse x. D’après une formule classique en 
théorie de l’élasticité donnant le cisaillement dans 
une direction quelconque, on a 

" K cos 0 


7 
v 


(cos20 — sin26). (15) 


Considérons d’abord le: cisaillement antagoniste 
créé par les accrochages directs accumulés au voisi- 
nage du point B. Si e est la distance moyenne de 
deux accrochages directs, le cisaillement dû aux 
accrochages directs est 


En 
= f Eu (cos28 — sin28), (16) 


Eq T° 


en désignant par 0, l'angle polaire. correspondant à 
l’accrochage’situé : l’extrémité du domaine d’accu- 


mulation. Il est évident que r et 0 sont liés par la 
condition 


S 


r sin(8—rx)= à ou r = — > (17) 


d'autre part, on a 


dx ÿ — dr 2 22 6 2 LAURE 
(dx =(dr}+r2(d8), dx sq dd oo 
D'où l’expression de (r:)£ 
K 3 sin? 2 0 sin 40 
NT 4 1 1 4 
Ge = ls 2 ne SU |. (19) 


Pour les accrochages inverses, on aurait de même 


: HÉAS sin? 2 0 sin 4 0! 
ANT Fr L ER Ai 1 É 
( Le PRES RS te | 

Si les domaines d’accumulation sont suffisamment 
étendus par rapport à Ô (et nous considérerons qu’il 
en est toujours ainsi), on a 0 et 07, 


(= 


d’où 
K 0 
Gh=rs(s-E): (20) 
Si nous admettons, pour simplifier, que les domaines 


occupés par les accrochages directs et inverses sont 
les mêmes (Na eu= N;ci= L,), nous aurons alors 


LES 


GNT (24) 


D'où, finalement, nous avons entre le nombre 
d’accrochages et le cisaillement antagoniste la rela- 
tion suivante : 


CSA r 
= (Th + (ri )e = A(SS RE Se =) (Na Ni). 


Cette équation est de la forme 7; = 6 (Ni — N) 
qui avait été admise par Kochendorfer et que nous 
avons utilisée systématiquement dans la théorie de 
la plasticité cristalline. 


Manuscrit reçu le 20 janvier 1948. 
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EFFICACITÉ DES COMPTEURS A RAYONS 
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a 


Sommaire. — Établissement de la courbe d'efficacité des compteurs à rayons y, à paroi d'Al, Cu, Au, 
dans le domaine d’énergie 0,3-3 MeV, par la méthode des coïncidences et des mesures absolues annexes. 
Utilisation des divers rayonnements des nuclides : 5V, AI, 2#%8Au, ‘Cu, Th C”. La connaissance de l’effi- 
cacité peut donner une présomption de la validité d’un schéma de désintégration encore imparfai- 


tement établi. 


Introduction. — [La méthode des coïncidences 
est fréquemment utilisée pour établir le schéma des 
niveaux d’un nuclide et la connaissance des effica- 
cités des compteurs pour les divers rayonnements 
permet de vérifier telle ou telle hypothèse d'un 
schéma inconnu ou de choisir entre divers schémas 
proposés. Le but de ce travail est la détermination 
expérimentale de la courbe d'efficacité des compteurs 
à rayons y, en fonction de la nature des parois et de 
l’énergie des rayons y dans un domaine allant 
de 0,3 à 3 MeV. 


Nous définirons tout d’abord l'efficacité d’un 
compteur à rayons y. Si un rayon y traverse le 
compteur, nous appellerons efficacité et désignerons 
par €, la probabilité qu’a ce quantum de donner 
lieu à une décharge. Soit une source radioactive 
émettant N rayons monochromatiques par unité 
de temps et donnant n coups par unité de temps dans 
le compteur considéré, nous appellerons rendement 
et désignerons par p, le rapport 


Nous pouvons relier le rendement à l'efficacité 
par la formule 


ie 
nl 
(w3 
© 


OV E+ 


où w est l’angle solide effectif et k le facteur d’absorp- 
tion dans les parois du compteur; ce facteur est de 
la forme 


3 : À à S SE. 
— est le coefficient d'absorption massique et ml'épais- 
F 

seur moyenne traversée en g : CM?. 

Le rendement est obtenu à partir de mesures en 
coïncidences et les formules employées seront déve- 
loppées dans le cas de chaque corps; w a été déter- 
miné à partir du rayonnement du Ra E, à l’aide de 
compteurs à rayons $, en extrapolant les courbes 


d'absorption dans l’Al obtenues avec un compteur 


\ 


à paroi de 1/1009 mm d'Al, un compteur à fenêtre 
de 194 d’AI et un compteur sans paroi (fig. 1). 


D DE D NE CR DEN EN IDR M 


us 


= 
ee 
æ 


Paroi compteur standard 


S 


nd 


Log activité unités arbitraires 


Epaisseur en u AL 
0 50 100 150 200 


Fig. 1. — Absorption dans l’aluminium du rayonnement. 
d’une source de Ra D + Ra E + Ra F en id utili- 
sation du compteur sans paroi. » | 


La masse superficielle moyenne traversée m a été 
déterminée par le rayonnement de dématériali- 
sation du nuclide #Cu et le rayonnement de 2,6 MeV. 
du Th C”, en utilisant divers compteurs à coques | 
de Cu, d'épaisseur variant entre o,1 et 2 mm: M 
m a été trouvé égal à la masse superficielle de la 
paroi du compteur. | 


Appareillage. — Compleurs, — Nous avons 
utilisé des compteurs CyHndriques à coque d’Al, 
de Cu et d’Au ayant 4o mm de longueur de fil a 
un diamètre intérieur de 17 mm (sauf le compteur ” 
à coque de Cu de 2 mm d'épaisseur et celui d’or É 
de 0,5 mm d’épaisseur qui ont un diamètre intérieur 
de 16mm); ïls sont remplis avec un mélange 


de ro.pour 100 d’alcool et go pour 100 d’argon, 
jusqu’à une pression totale de rocm de Hg. Les 
coques sont, soit en Al : épaisseur o,r et 1 mm, 
soit en Cu : épaisseur 0,1, 0,2, 0,5, 1, 1,5, 2 mm, 
soit en Au : feuilles d’Au, d'épaisseur 5 ps, 10 pet 
6,5 mm enroulées dans un compteur d’Al-de o,1mm 
d'épaisseur. 

. Pour la détermination de l'angle solide effectif, 
par extrapolation de la courbe d’absorption du Ra E, 
nous avons utilisé un compteur en laiton de mêmes 
dimensions portant une fenêtre de 1 cm? fermée par 
une feuille d'Al de 19 p. Le compteur sans paroi, 
utilisé pour la même mesure, est un compteur forme 
cloche, non fermé, de 20 mm de diamètre, l'extrémité 
libre du fil est située à 20 mm de la source, l’ensemble 
est placé dans un dessiccateur rempli du mélange 
utilisé pour les compteurs précédents. Un disque 
porte-écrans intercepte le rayonnement émis de la 
source vers le compteur, le disque .porte sur sa 
circonférence trois masselottes de fer doux qui 
permettent de commander sa rotation extérieu- 
tement par un aimant. Le disque porte six écrans 
ayant comme épaisseur respectivement : 218, 104, 
58, 21, 6 et op d’AI (fig. 2). 


Fig. 2. — Compteur sans paroi. 
a, compteur; b, disque porte-écrans; €, masselottes de fer 
- doux; d, robinet permettant le pompage et le remplissage 
du compteur; e, source. 


Supports de compteurs. — Les compteurs peuvent 
être fixés en différents logements situés à des dis- 
tances fixes de la source, permettant diverses combi- 
naisons tout en évitant des erreurs de repérage des 
distances (fig. 3). 


Sources. — Les sources sont soit déposées sur une 
feuille d’acétate de cellulose (épaisseur + 8), 
tendue sur un cadre de plexiglass, soit, dans le 
cas de 28A1 et 18Au, collées directement sur le cadre 
de plexiglass par une trace d’une solution d’acétate 
de cellulose dans l’acétone. 


Appareils de numération. — Les compteurs sont 


EE 
| 


6 EFFICACITÉ DES COMPTEURS À RAYONS y 43 


reliés à un amplificateur à quatre coïncidences, 
à une échelle de 1000 et à une échelle de r00 cons- 
truits par le Commissariat à l'Énergie atomique. 


Tension fl compteur 


cu 


Ù Fig. 3.— Support compteurs. 

S1, So, S3, emplacements possibles des sources; A, B, C, D, E 
E, G, H, I, J, emplacements possibles des compteurs; 
e, e’, e”, etc., emplacements possibles des écrans. 


Ces appareils permettent la mesure simultanée du 
nombre d’impulsions dans chaque compteur et du 
nombre de coïncidences. Le pouvoir de résolution 
est de l’ordre de 0,5.10o7fs en coïncidences doubles, 
pour toutes les combinaisons de compteurs deux 
à deux. Le taux de 10000 impulsions : mn dans 
un compteur n’a pas été dépassé (sauf dans les 
mesures sur Æ8Au) dans le but d'éviter les pertes 
dues au temps mort du compteur. De ce fait, les 
coïncidences fortuites sont négligeables (sauf dans 
les mesures sur SAu où elles ont été évaluées). 
Le bon fonctionnement des diverses combinaisons 
de compteurs en coïncidences a été testé par la 
méthode décrite par Norling [r]. Ce test a été effectué 
à diverses reprises au cours des mesures qui ont, 
de plus, été encadrées par des mesures d’étalons. 


Mesures. — Nous avons utilisé des corps dont 
le schéma de désintégration a déjà été étudié : 


19 SV, sous forme de V,0;, pour l’énergie 
de 1,46 MeV; 


20 28A], pour l'énergie de 1,8 MeV; 


30 24Na, sous forme de CO,Na,, pour les énergies 
de 1,38 et 2,76 MeV; 
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4° M8Au, pour les énergies de 0,3 et 0,4 MeV; 
. 50 61Cu, sous forme de CuS, pour l'énergie du 
rayonnement d’annihilation de 0,51 MeV. 

D'autre part, des mesures directes ont été effec- 
tuées avec : s 


19 Th C”, d’une source de R Th, pour l’énergie 
de 2,6 MeV; 


20 61Çu, pour le rayonnement de 0,51 MeV. 


Les résultats des mesures sont donnés dans le 
Tableau I. 


TABLEAU I. 
k Efficacité (>< 104). 
Energie a  — 
Nuclide. (en MeV). Cu. Al. Au. 
CN RE A0 LOMME EUR D PTE 12/07 
SAINTES MID 00,7 — 192 0220.09 
Total 27,0 0,5 198;8E0,4,: 29/22 
2ANa2 09-70 19 18,8 LOT 
1,38 8,5 NON ne 
Total 4,700, 1 NERO AD DO 
198AU. ; {0,3 2 IA He 
0,4 2,70 2,9 
é Cu o,510b JD =ONL 9, J==0 2 SI OMEET 
LH OLD rOS ML OS 0.0 0,0 74 
1bCT 25% ,;02 LOS ONDES, 4 0,0 RTE 0,D 
“Valeur moyenne. — ?Mesure en coïncidences, — ‘Mesure 
directe. — dValeur estimée. 
Vanadium. Nous avons utilisé un porteur 


de 2 mg de V;,0, en poudre, irradié pendant ro mn 
par neutrons lents. Le schéma de désintégration du 
nuclide 52V formé a été étudié par T. Yuasa [2], 
J. Martelly [3] et par R. Bouchez et nous-mêmes [4]; 
on peut le représenter par le schéma 4 a. Dans «la 


pris ass 
p 33 FY04 


y/8 Y28 


28n1 2 198, 
b C 


Fig. 4. — à, b, c, d. 
a, SV; b, #Al; c, 2Na d, Au. 


mesure des coïncidences B—#7, en utilisant un 
compteur à rayons 8 dont la paroi arrête les élec- 
trons de la transition isomérique (3/10 mm d’Al) 
les négatons du spectre B seront seuls comptés. 
La proportion de rayons y donnant des décharges 
dans le compteur à rayons $ est de l’ordre 
de 1 pour 100 et peut être négligée à la précision 
des mesures; d’autre part, l'épaisseur du compteur 
à rayons y est suffisante pour arrêter tous les élec- 


trons. On peut écrire. 


Ns =Nes, 
Ny = Ne), 
Ngy= Nopey, 


N8, N,, N8_., N, étant respectivement le nombre 
d’impulsions provenant du compteur à rayons f, 
du compteur à rayons y, en coïncidences et le . 
nombre de désintégrations de la source pendant » 
l’unité de temps; pg, p, les rendements des comp- # 
teurs 8 et y. Nous avons, avec les notations de . 


l'introduction, 
ÊYES wok, 
d'où nous pouvons tirer la valeur de l'efficacité. 


Aluminium. — Le porteur est une feuille d’Al 
de 3 cm sur 1 cm, de 4 x1o-? mm d'épaisseur, pur 
à 0,997, irradiée pendant 10 mn par neutrons lents. . 
Le schéma de désintégration du nuclide Al peut 
être représenté par le schéma 4 b [5] et le rendement 
et l'efficacité sont obtenus en utilisant les mêmes "| 
formules que dans le cas de 52V. 


Sodium.e — Nous avons irradié du CO;Na, par . 
des deutons de 6,8 MeV pendant des temps allant 
de r à 5h. Le produit est dissous et neutralisé. 
et une goutte de la solution déposée sur un porteur « 
d’acétate de cellulose, constitue la source après 
séchage à l’étuve (). La source décroît avec une. 
période pure de 14,8 h à 1 pour 100 près. En admet- 
tant le schéma proposé par Elliot et Alburger [6-8], 


schéma 4 ce, et en négligeant les électrons de conver- 


sion, nous obtenons une valeur moyenne du ren- 
dement, d’où, le point obtenu pour le V servant 
de comparaison, nous pouvons tirer les valeurs + È 
l'efficacité pour les énergies de 1,38 et 2,76 MeV. 
Nous avons utilisé pour déterminer le mel 
les formules (avec les notations habituelles }) : 


N8 — Ne8, 
Ny = Noy(skæ+ 4) = Npyiotat, 
Nyy= Nogwy(sk+ek), 


:, €, k, K', étant respectivement les valeurs de 


l'efficacité et du facteur de transmission pour les 4 


raies de 1,38 et 2,76 MeV. 


Or. — Nous avons utilisé différentes feuilles d'or - 
de 3 cm sur r em et de 1 cm sur 1 cm d’épaisseurs 10, M 
5 et 0,5 p, irradiées par neutrons lents ou par deutons : “ 
pendant des temps variant de quelques heures à 
un mois. Nous avons trouvé, dans tous les cas, 
une période pure de 66 h à 1 pour 100 près. Le schéma M 
le plus probable d’après nos mesures en coïnci- 
dences [17] est celui représenté par le schéma 4 d. 


(:) Nous remercions M. P. Süe et ses collaborateurs de la 
préparation des sources de Na. 


(4 |En négligeant les éctrens de conversion et la possi- 
_bilité d’une raie y faible de 2,5 MeV, on obtient 

lle rendement en utilisant les mêmes formules que 
. dans le cas de 24Na. 


Cuivre. _ Nous avons utilisé la réaction Ni (d, p) Cu 
par la méthode décrite par Savel [o], pour obtenir 
- le nuclide S!Cu. Après séparation chimique et centri- 
: fugation, le précipité de CuS est déposé sur un 
| d'acétate de cellulose. L'activité du produit 
_ mesurée quelques heures après la fin de l'irradia- 
pe décroft avec une période pure de 3,4h à 

5 pour 100 près. La source est entourée d’une 
| fouille de Pb de 1 mm d’épaisseur; les B+ sont ainsi 
arrêtés. En faisant l'hypothèse, que deux rayons y 
! de 510 keV sont émis à chaque annihilation, et qu'ils 

_ sont émis dans deux directions diamétralement 
| opposées, les mesures en coïncidences ÿ—7Y per- 

mettent de remonter au nombre de désintégrations 
… et à l’efficacité. On a les formules 


Ny —= 2Nwe,k, 
N}, — 2Nwe, ko, 
Ny-,= Nw s$ki ko, 


_N,N,,N,, N;_;, étant respectivement les nombres : 
- de désintégrations, de décharges dans le compteur 1 
. et le compteur 2 et en coïncidences; w l’angle solide 

effectif des deux compteurs (le même pour les deux 
compteurs); les autres notations sont conformes à 
- celles déjà mentionnées dans l'introduction. Remar- 

quons que la corrélation angulaire des deux y donne 
des formules différentes de celles établies dans le cas 
général des coïncidences y — +. 


Mesures absolues directes. — Nous voulons 
parler des mesures absolues que nous avons effec- 
tuées à titre de contrôle avec le rayonnement 

dun Th CG'et.le rayonnement d’annihilation du 
_ nuclide 6!Cu. 


Thorium C". — Nous avons utilisé le rayonne- 

. ment, filtré par 3,5 cm de Pb, d’une source de R Th 

équivalente à 0,126 mg de Ra préparée et étalonnée 

par la méthode du grand condensateur à l’Institut 
du Radium (?). 


Cuivre. — Le nuclide $Cu a été préparé par la 
- même méthode que précédemment, nous nous 
* sommes servis du même échantillon que celui 
: employé dans les mesures en coïncidences ; nous avons 
utilisé un compteur à rayons B en Al de :/r0° de 
millimètre d'épaisseur et le compteur sans paroi 
précédemment décrit, pour déterminer l’absorption 
dans la paroi et en utilisant un diaphragme cana- 
lisant un faisceau de B*, nous avons pu déterminer 
le nombre de désintégrations de la source. De là, 


(:) Nous remercions M. G. Bouissières et M° Chamié 
pour la préparation et l’étalonnage de cette source, 


. (Tableau 
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on obtient directement le rendement et l'efficacité. 
Cette méthode est peu précise, ce qui peut expliquer 
les différences avec les résultats obtenus par la 
méthode des coïncidences. D’autre part, dans les 
mesures en coïncidences, nous avons supposé que le 
rayonnement d’annihilation n’est constitué que de 
quanta de 0,51 MeV, alors que la théorie prévoit un 
rayonnement de 1 MeV (de l’ordre de 10 pour 100 du 
rayonnement total dans le cas des éléments lourds), 
et que ce rayonnement a déjà été observé [ro]. Des 
expériences sont en projet à ce sujet. 


Résultats. — Nous avons comparé nos résultats 
I) à ceux obtenus par Norling [r], 
Bradt et al. [11], Maiïer-Leibnitz [r2] et Dun- 
worth [13], ainsi qu'aux résultats théoriques et 
semi-théoriques de Von Droste [14], Yukava et 
Sakata [15] et N. Marty [16]. Les courbes 


sont dessinées (fig. 5, 6 et 7) pour le Cu, l’AI et l’Au. 


Fig, 5. — Efficacité des compteurs à coque d’Aluminium, 
Maier-Leibnitz, + 2220 —; Bradt et al, ------ 
N>Marty,- 14 ; Renard, o, mesures en Caïncidences: 


e, mesures absolues directes. 


Le décalement systématique de certaines courbes 
provient vraisemblablement de différences entre les 
dispositifs expérimentaux. La figure 8 reproduit 
nos résultats quant aux efficacités pour les divers 
matériaux utilisés. La figure o représente les rap- 
“ en fonction de l'énergie du rayon- 


2 
F 


ck 
ports 7 — 


NL ARE 


Fig. 6. — Eficacité des compteurs à coque de cuivre. 
Maier-Leibnitz, — + —. — Brodt bal ee ; N. Marty, 


e, mesures absolues directes. 


s 
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EYy en MeV 


Fig. 7.— Efficacité des compteurs à coques de Pb et d’Au. 


--, Au: Renard, 
© mesures en Coïncidences; e, mesures absolues directes. 


Pb : Maier-Leibnitz, — : — + — + ; Bradt et al, 


nement; les facteurs ket k' sont calculés pour une 
épaisseur de paroi égale au parcours maximum des 
électrons secondaires pouvant être créés dans 
cette paroi; nous avons utilisé la valeur moyenne 
déduite des relations parcours-énergie des rayons B, 


de Coryell () et de Flammersfeld [r8], ainsi 


(°) E, SR 


causerie). 


R + 0,094 
0,526 


(formule donnée au cours d’une 


| Ey 1030 


que d une ile provenant de Princeton (9. 


Ces courbes ont été établies en vue de vérifier les 


| efcncités obtenues par la méthode des coïncidences, 


02 03 04 05 06 08 1 
Fig. 8. — Efficacité des compteurs à rayons Y. 
mr ra , Coque d’Al; 
, Coque de Cu; 
D COQUE AU. 


dans l’étude d’un nuclide dont le schéma de désin- 
tégration n’est pas encore définitivement établi ou | 
est controversé. Ces courbes se sont révélées utiles 
pour notre étude du schéma de désintégration 
de 1*#Au [17]. Dans cette étude, nous avons supposé 
la continuité de l'efficacité avec l’énergie; en fait, 


( Exemplaire fourni par M. Magnan. 
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A 


dans la région utilisée, l'éventualité d’une discon- 
tinuité est improbable théoriquement, ce qui justifie 
notre manière de procéder. 


À 
: 
= 
: 


nl 
Fe 
His 
FF 
Fi 
ie 


2 
AL 
IL 


ot 
UE 


0.8 = 
de ere. 
"O2 03 04 05 6 08 1 CRE 7 
Ey en MeV 
Fig. 9. — Rapports des rendements des divers compteurs 


en fonction de l'énergie des rayons y; pour une épaisseur 
de paroi correspondant au parcours maximum des élec- 
trons secondaires. 


Nous tenons à remercier M. le Professeur F. Joliot, 


.qui a dirigé notre travail et nous a aidé de précieux 


conseils. Nous exprimons notre gratitude au C.N.RS. 
pour l'allocation qu'il nous a fournie. 


Manuscrit reçu le 21 mai 1948. 
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SUR QUELQUES PROPRIÉTÉS DES CORPUSCULES DE SPIN - 
DANS LES CHAMPS ÉLECTROMAGNÉTIQUES 


Par GÉ£RrARD PETIAU. 
Institut Henri-Poincaré. 
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Sommaire, — L'introduction dans l'équation d'ondes de l’électron de Dirac, représentant d’une 


Re UE : En ee 
façon générale, les corpuscules de spin Her Le d'un terme d'interaction électromagnétique s’exerçant 


directement par les champs et non plus seulement par les potentiels, permet de représenter un corpuscule 
possédant un moment magnétique propre non lié à l'existence de la charge. 

Nous examinons ici quelques propriétés de ce corpuscule que l’on peut déduire de Pénéation d'ondes. 
Après avoir défini un opérateur vitesse généralisé, nous montrons comment l'introduction d’un tenseur 
d’impulsion généralisé permet d'étudier les forces résultant de l’action du champ sur le corpuscule, 
celles-ci se décomposant en une force correspondant à la force de Lorentz classique et en une force 
due au moment électromagnétique propre ou force de spin dont nous discutons la possibilité de mise 


en évidence. 


Les résultats obtenus sont ensuite retrouvés dans l'étude de l'approximation non relativiste de 
! l'équation d’ondes considérée, d’abord dans le cas du champ électromagnétique quelconque puis, plus 


en détail, dans le cas du champ électrique central. 


1. La théorie de l’électron relativiste de Dirac 

; 5 Ne 
représente l’état du corpuscule de spin 2 on 
l’absence de champ extérieur par les fonctions 
d'ondes dx (k = 1, 2, 3, 4), solutions du système 


d'équations linéaires aux dérivées partielles du 
premier ordre : 


[Po + p.x + Mmocai ]bx = 0, (1) 
où 
0 ro) É h 
po= it TL pre h —= + 


x, «, représentant un système de quatre matrices 
du quatrième rang, telles que 


du y ay = 2 dy CRT ESS En 0e 
Lorsque ce corpuscule portant une charge € se 
trouve placé dans un champ électromagnétique 
caractérisé par les champs électrique et magné- 
tique E, H et par les potentiels scalaire et vec- 
teur À,, A, l'équation d’ondes déterminant les 
fonctions Ÿ, se déduit de (1) en remplaçant les 


opérateurs ps, P Par Po + = As p + = A, ce qui 
revient à compléter l’équation (1) par l'expression 
invariante dans une transformation de Lorentz 


€ 
el4+ A.x| ( 
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et à écrire celle-ci sous la forme 
Po + p.x+ : [Ao+ A.a]+ 70 ca | VE = 0: (3) 
L'interprétation de cette équation conduit à 


attribuer au corpuscule un moment cinétique 
propre ou spin et un moment électromagnétique 


propre auquel correspond l'énergie potentielle 


A LE. .r + Hp] (4) 
en désignant par x et pm les matrices vecteurs de 
composantes 


Rp= Tapas, Up= Spa = 14 ar, 


CP; JET; 2; 3). 


Cependant, W. Pauli a montré (1) que l’expression 
invariante (4) pouvait également être introduite 
directement dans l’équation d’ondes (1) qui s’écrit 
alors 


Sp = lag ar 


[po+ pa + [Ex + H.p] + moca|dr= Cut) 


et dont les solutions représentent les états d’un 
corpuscule non chargé mais possédant un moment 
magnétique propre {4. 

Cette équation a été considérée depuis par de 
nombreux auteurs pour représenter le proton ou le 
neutron selon que l’on tient compte, simultanément 
ou non, des deux termes d'interaction (2) et (4). 
En particulier, cette équation a été utilisée dans 
des travaux récents de Caldirola (?), G. Breit (5), 
R. G. Sachs (#). Signalons, en outre, que H. Bethe (5) 
s’en est servi pour représenter le neutrino dans 
l'hypothèse où ce corpuscule non chargé, de masse 


@) W. Pauzr, Handbuch der Physik, 24/4, p. 533. 

(2) GALDIROLA, Relativistic correction in calculating ‘the 
magnetic moment of deuteron. Phys. Rev., 1946, 69, p. 608. 

@) G. Brerr, Relativistic corrections to magnetic moments 
of nuclear particles. Phys. Rev., 1947, 71, p. 100-402. 

(*) R. G. Sacus, On the magnetic moment of deuteron. 
Phys. Rev., 1947, 72, p. 91-98. 

(5) H. BETHE, Tonisation power of a neutrino with moment 
magnetic. Proceed. Cambridge Phil. Soc., 1935, 84, p. 108-115, 


très petite, 


NE RE ut DIET 
. BOPESRET er DR 
AL AE sEAR) fers 
; 


ionisation des atomes par l’action directe d’un 


. moment électromagnétique propre sur les électrons 


intraatomiques. 


En présence de ces travaux, il nous a semblé 
qu'il serait utile de préciser quelques-unes des pro- 


_ priétés que l’on peut déduire directement de l’équa- 
_ tion d'ondes du corpuscule considéré. Nous nous 


proposons d'exposer ici les résultats que nous avons 


. obtenus. 


2. Nous considérons le corpuscule dont les fonc- 


tions d’ondes Ÿx (que nous écrirons simplement d 


É à l’avenir) sont solutions dans le champ électro- 


magnétique À,, À, E, H de l'équation d’ondes 
Lpo+ pa+g[4o+A.a] 


— ÎLE :x + Hp] + mocai|d = 0, (6) 
g et j étant deux constantes caractérisant les inter- 
actions. 
A l'équation (6), nous associons l'équation 
conjuguée 
Po" +p"4"a 
4" “e(40+ Ac) 


_—f(E.x + H.p) + moca] = 0. (79) 


Multipliant ces équations par «,, nous en déduisons 


la relation 
moe | y* 4 de 
+ f #Lpiai+ pl de 
—n [ ST sud] de 


3 JT 
e f v*vae _ 


mc fyta de 
7 h 0 


+2 f va do—2f f #{H.sl dé = 0. 


. Tenant compte de la condition de nullité des Ÿ 
-aux limites du domaine d’intégration et de l’hermi- 
ticité des opérateurs, cette relation se réduit à 


Jto+ g A6) — f(H.0)]4 de 
(8) 


en posant 
(9) 


De même, en multipliant (6) et (7) par ps, 
nous obtenons, par combinaison, la relation intégrale 


Lo = (po + 8 Ao)xi — f(H:5). 


2 fuf@+sma— le A s)— He] | de = 0 


(o = iagar, ldrap, ap@g, O0 = ap Xg@r). 
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serait susceptible de provoquer une 


(10) 
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Or, nous avons, en l’absence de champ extérieur, 
Ë I 
Udg = 


_ * * { 

e [y Je Por: Ÿ de, 
d ñ 
‘22 AS EAN ARTS ANSE 
SJ? [x D a]va - fo ne 


D'autre part, quel que soit le type d'interaction 
considérée, nous avons 


(11) 


dE 3 
A Yixvde=— cf ÿ'aÿ de. (12) 


dre 
Combinant cette relation avec (10), nous en dédui- 
sons 


dise le 
——- DT RE 7 | d de 
dé : 2 MoC : 


I 


m0 


q* Lo + gA)o: 
Le f[TE \6]l— H || Ÿ de 


ï . 
=— — d* IL do 
1110 ee ‘ 


(413) 
avec 


N=(p+gA)u—/[[EAs]—-Hal (14) 


Les relations (8) et (13) généralisent (11), le vec- 


teur Il, Il correspondant au vecteur p, pv —0p P. 
0 Ù m 


de la Mécanique ondulatoire sans spin. 


3. Nous avons déjà rencontré le quadrivecteur IL,, 
Il dans un travail précédent (f) où nous avons 
étudié comme cas particulier les équations d’ondes 
du second ordre correspondant à l'équation (6). 
Par itération de celle-ci, on obtient facilement 
((po+g 40) —(p+gA)— mic 
— 2/(H.h)(pe+ 8 40) 
+2f([E A8]l.(p+sA) 
+2/(H{(p+gA))ia;—/f'(E— H:) 
— 2/2 E.H)io+ À g[(E.ix) —(H.s)]}Ÿ = 0. 
Cette équation, dans laquelle nous avons posé 
J As = AjyXo ds TX 
n’est pas hermitienne, mais par multiplication 
par «,, elle devient hermitienne et s’écrit 
a [(po+ #40) 7 (pP+sA)— mic] 
— 2f(H.5)(po+ 8 40) 
+ 2f([E À s].(p+#A)) 
— 2 f(H(p +gA)) So 
— PTE Ha 
— 2f2(E.H a; — Ag[E.x + H.p| 1% — 0: (M9) 


Introduisant dans cette équation les expressions (9) 
et (14) de I, et Il, nous obtenons, après division 


(6) G. PETIAU, J. dé Physique, 1947, 8, p. 134-139. 
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Uo— Mail — mc; ] 


J” [CE — H° co; + 2(E.H )c; | 


Cette équation met en évidence, d’une part 
l'interaction électromagnétique due à la charge 


. électrique et au spin 


hs | 
Smet ul À 


et, d'autre part, deux termes d'énergie potentielle 
correspondant à l'interaction en f, l’un classique : 


f° 


270 


[E?— H° Je, 


l’autre 
nue (E.H)c; 


mo 
purement quantique et caractéristique des mélanges 
d'états à énergies positives et négatives. 
4. L'équation (6) écrite sous la forme 
Ld=o (17) 


nous permet de calculer la dérivée par rapport au 


temps des grandeurs caractéristiques du corpuscule 


au moyen de la formule de dérivation quantique 


d * ic *r 1 A) É 
si AŸ de = = [4 [ZA AL] ar, (18) 


A étant l’opérateur correspondant à la grandeur 
considérée. 


En particulier, si nous considérons l’opéra- 
teur p + gA ayant la variance de la composante æ{ 
d’un tenseur du second ordre, nous obtenons faci- 


lement la relation 


d à 
n J Y@+sa)vas 
= ge [ YIE+HA «jy de 
—fe [4 grad(E.z-+H.p)bde. (19) 


Le premier terme du second membre peut se 
transformer en utilisant les équations (6) et (7) 
multipliées par 

Ec;+{H / «lo 
et par 
Ec, + ci | H /\ œ |. 


On obtient alors, par combinaison, la relation 


intégrale 
ee À a]u de 


= 2 fr ilEur+ #40) 
—[HA(p+sA)]s 
é — fe [4 de 
J'Y A 18 46 


fy Brad(E SE AR Ho 5 és Gus : 


h d 
2mM0C? dt 


2mM0C 


si nous tenons compte de cette A SE dans (9), 
celle-ci devient 


d * 
a? 
= À qe Has 
—[HA(p+gA)]« 
à ë USE: a 


+[ fe] f grad(E. x + H.p)y de. @1) 


p+sA 1]uae à 


Mais, d'autre part, les définitions (9) et (14) de I, * 
et IL et l’identité vectorielle :. +3 


[A AB A CI] = (B(A.c) —(C(A.B) 
nous donnent - à à 
Ell—[H A1] | 
= [E(po+ 8 40) 


[HA (p + #A)i]a JE. H)s.*_ (22) "4 
Par suite, la relation (21) s'écrit encore 4 
d 
nf Y{p+sa+ e EUR À») | v de 
= E [4 [EL TH AM] y de 


É Es — je) [ d*grad(E.z+ H.u)bdo (93) 


Cette relation généralise celle que l’on obtient (°) 
dans le cas de l’équation de Dirac usuelle où l’on a 

£= =) f=0 
et qui s'écrit 


Fe 
a Y [p+ia c à /\ =1]4 er 
Er” sal We [E(ro+ g40) | 
—[HA(p+gA)l|uvde 


ñ 
ce [+ grad(E.r+Hip)ÿde (24) 
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(7) G. PerrAu, C. R. Acad. Sc., Paris, 1948, 226, p. 313-314. 


mettant en évidence au second membre la force 
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de Lorentz classique et la force électromagnétique | 


due à l'existence du. spin ou force de spin corres- 
_ pondant à l'énergie électromagnétique propre 


= : 
eh 


He [E.x + H.u|. 


Néanmoins, l’équation (23) qui met en évidence 
au second membre la force de Lorentz généralisée 


_ représentée par l'opérateur 


_ avec, pour coefficient d’une façon ‘additive, les 


moments 


À (Elo [H An 


_et la force de spin d’opérateur 


erad(E.x + H.y) 


en 


_ satisfaisante. En effet, l'opérateur p + gA peut 


être complété par des opérateurs de même variance 


_ construits au moyen des grandeurs du champ et 


des matrices caractéristiques du corpuscule. En 


particulier, nous pouvons compléter cet opérateur 


par l'expression —f(E /w) en considérant pour 
tenseur d’impulsion totale l'expression 


p+sA—/J[EAR] (25) 


Afin d'examiner les modifications que le choix de 
_ cette expression entraîne dans la représentation des 
. forces appliquées au corpuseule, nous devons dériver 


_ l'opérateur [E À pe]. 


En utilisant les équations (6) et (7) on obtient 


facilement la relation intégrale 


en J'TE AIYde 


c dt 
= frfate ne srd(e.s |y de 
+? [AE A (+44 
[HA s])| 4 dy. 


Remarquant que nous avons 


[an 


(26) devient 


Éd: ; 
L* d de 
e a FE A us su 


jy grad(E.x + H.p.)v de 


—=— grad(H.p) + (pgradH), 


+ [ y (.gradH)v de 
& [uw [E À ((p+ gA)s 
+ /TH À s])| ly de. 


Alors que, dans (21), nous avions introduit 
N=(p+gA)x—/f([E As] —Ho), 


(28) 


et —fc, n’est pas complètement 


to 


nous sommes conduits ici à introduire le vecteur dual 


N'=(p+gA)a;+ /f([H A 5]+Ec). (29) 
De même, nous introduisons le dual de 
Ho = (po+ gAo)ai— f(H.0) 
par 
6 = (po+ 8 Ao)a; + f(E.0). (30) 


Or, par combinaison de (6) et (7) multipliées par H «:;, 
nous obtenons la relation intégrale 


à Je Hi + gAna+ fo] de 


= f v(wgrad Hyde (31) 
et, par suite, (28) devient 
fa f VIE Aa de 
= 2% [yen -[e Ant] 
re [y grad(E.x + H.y)Ù de. (32). 


Complétant (23) avec cette expression, nous 
obtenons finalement la dérivée de l’impulsion géné- 


- ralisée (25). : 


: | 
n fyp+-sa-reau 


Len L de 
+ els /\ s]]v de 


ee 


[Em —[H AN], de 


2 fc . À à 
+2 fo[Hnt—[e À mi]va 
Ji 

gen 
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fe Srd(E 2 Hidden (os 


Nous voyons apparaître au second membre deux 


- « forces de Lorentz » duales l’une de l’autre, asso- 


on x € x 
ciées, l’une à la charge g=—<) l’autre à une 


charge fictive s' que l’on peut déduire de la cons- 


tante f par la relation &' — TE 
force de spin est alors celle qui correspond à la 
charge vraie, c’est-à-dire à la constante 9, résultat 
en accord avec celui que nous déduisons de l’équa- 
tion d’ondes du second ordre (16) où le terme 
d'énergie potentielle dû au moment magnétique 
propre ne dépend que de la constante g. 


Remarquons que si nous posons g° — 2 le ten- 


f. Mais la seule 


_seur d’impulsion généralisé intervenant au premier 


membre de (33), s'écrit encore 


o #, F o| — 2" fs / al ; ( 34 
PS A me lélH A #1] [E À will Ne. 
Dans le cas où le champ extérieur est uniquement 


un champ magnétique indépéndant du temps, les 


999 
formules (33) ou (24) se réduisent à 
: ol y [p+s 
SE Le H À (p+ gA)a]$ de 
+ [Se =#| [+ grad(H.p)d de. (35) 


[H A s]]ya 
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Dans ce cas, à côté de la force de Lorentz pro- 
prement dite, les forces en g et f correspondant au 
moment magnétique propre interviennent d’une 
facon additive. 

Si l’on examine les conditions dans lesquelles on 
peut évaluer la force de Lorentz représentée par 
l'opérateur 


H /\(p+gA)x| 


on voit que cette évaluation exige la mesure simul- 
tanée du champ gH et de la vitesse (p + gA)=— “= 


Or, d’après la théorie générale de la mesure 6, 
on sait que si deux grandeurs physiques observables 
sont représentées par les opérateurs linéaires et 
hermitiens À et B, on a 

AA ABS - [A4 BT, Bye 
[A, B] représentant le commutateur de À et B 
et AA, AB les écarts quadratiques définis par 
(AAÿ=(4- 4),  (ABÿ=(B—B). 

Appliquant ce théorème, nous voyons que l’éva- 
luation de la force de Lorentz introduit une erreur 
supérieure à 


gh \ 
— | grad Ho, | ou; 
20 


Nous retrouvons ici un résultat qui a fait l’objet 
de nombreuses discussions (°) d’après lesquelles on 
conclut généralement que la force de spin ne peut 
être mise en évidence par une expérience directe, 
l'incertitude dans la mesure de la force de Lorentz 
étant toujours supérieure à la force de spin (ou de 
moment électromagnétique propre). 

Mais nous voyons ici que ce raisonnement ne 
porte que sur la force de spin dépendant du facteur g 


( c'est-à-dire de la charge si g — :) + Si le coefficient 


d'interaction électromagnétique directe fe est nota- 

blement supérieur au moment magnétique propre 
k ch : 

_ = » la force de spin pourra 
mo 2.0 C 


être mise en évidence par une expérience directe. 


(5) Voir, par exemple : L. DE BroGLre, Théorie générale des 
particules à spin, p. 30-32. 

(*) Voir, par exemple, le Rapport de M. W. Pau au 
Congrès Solvay de 1930 sur le Magnétisme. 
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Nous retrouvons encore ici une des conclusions du ‘4 
rapport cité de M. W. Pauli (1). 

On sait, par ailleurs, que M. J. Thibaud (1), 
étudiant le comportement de certains rayonnements 8 
dans des champs magnétiques inhomogènes a été 
conduit à admettre l'existence de particules très 
pénétrantes caractérisées par un moment électro- … 


magnétique important et qu’il a désignées sous le 4 


nom d’électrinos. Si l’on admet que ces particules 
satisfont à l’équation d’ondes de Dirac, la considé- 
ration d’une interaction magnétique directe en } 
[analogue à celle introduite par H. Bethe (5) dans 
le cas du neutrino] donnant un moment magnétique 
prépondérant devant celui en g permettrait d’exa- 
miner indépendamment l’action de la force de 
Lorentz et celle de la force de spin correspondant au 
moment magnétique en f et de justifier, par un 
équilibre entre ces forces, les répartitions corpuscu- 
laires observées. 


5. Nous allons retrouver les résultats précédents 
en examinant l’approximation non relativiste corres- 
pondant à l'équation (6). 

Pour cela, suivant le procédé introduit par 
W. Pauli et repris par R. G. Sachs, nous allons 
utiliser la représentation donnée par Dirac pour les 


matrices « en écrivant 
x — GP1, a, = 3 


les s et les p étant deux systèmes de trois matrices 


à deux lignes et deux colonnes, indépendantes, 
tels que | a 
SpSq = — Ü9r PpPq=— {Pr . CpPq=PySp 


(P;-9, T =1,2, 3) 


et pour lesquels nous prenons la représentation de 
Pauli : 


c1 OST Go) (0) t 
où pi Ho: 43 ou pe EE) è 
O3 x O 
net : 


Nous avons alors 


tp Lg = Or, LA pQq Ur = C0 = Pi) Lp Ag Ar ay = y = P9, 


tapaq ui = Ur = OrP3) LapU = Tp = —Cppa 


Si, au lieu de repérer les quatre fonctions d’ondes de 
par l'indice k = 1, 2, 3, 4, nous introduisons deux 
indices { = 1; 2, m1, 2 écrivant W,,, l'indice / 
correspondant aux matrices a, 
matrices p, les fonctions d.,, VW, s’écriront d/, les 
fonctions d,, d, s’écriront 4. Nous écrirons, dans 


ce qui suit, ©, pour 1, 9 pour 7, en omettant 
l'indice L. 


En mettant en évidence ces fonctions Pr Et Do, 


(Gt) MY: Hu C.1R;: Acad. Sc. Paris, 1946, 223, 
P. 984; 1947, 224, p. 739 et 914. 


l'indice m aux * 


Er 


as Tr 


| LS 
| Ÿ “a We, 
ENG. 


__ QUELQU 
| l’équation (6) se décompose en donnant le système 
Ë [po+gAo+ moc— f (H.c)]o: 

+ ((p + gA + ifE).s)pr=— 0, 
[po+gAo— moc+ f(H.0)]9: 
+ (p+gA—if#).c)o= 0. 


| 


(36) 


Si, dans l'opérateur d'énergie cps, nous séparons 
l'énergie massique propre mo c? en écrivant 

| Po= Po + MC, 

cp, correspondra à l'énergie non relativiste. 

Si, en outre, nous avons 

8 Ao— f H < mc, 

W 

(6 


Po À 


avec WÆ mc? 


la première des équations (36) nous donnera 
| approximativement 


es ((p+gA+ifE].c)e: (37) 


20 C 
| tandis que la seconde des équations (36) s'écrit 
D  [po+g40+/f(H.s)le 
É +[p+gA—ifE#£).e]s1=0. (38) 
Reportant (37) dans (38), les fonctions 9, seront 
. déterminées par l'équation 


pig Au JUL) 


I 
2 Mo C 


X[P+sA+ifE).s]lsi—o. (39) 


[(p+gA—1/E#).0]] 


Æ Or, nous avons, tenant compte de div E = 0, 
(P +gA—ifE)(p+gsA+ifE)=(p+gA)+/f?'E? 
et en utilisant l'identité 

| ([A A B].C)=—([A A CI-B) 

| A [rp+ &Ap— if Ep][p1+ 8 44+ Ü Es] 

EP s 
— [Pat 8 Ag — if Es][pr+ 84» + if Ep]| SpSq 

= £h(H:6)+ 2 f (CE /\ 6].(p + gA)) + ifh(RotE. ), 


__ L’équation (39) devient alors 


2mMm0C 


pit g Ait J (He) — 


x [(p+sA} 
) +2. f([E Ask(p+ gA)) + f2E° 


+ gh(H.s)+ if h(RotE.c)| 


9»=0. (40) 


Posant s 
P=p+gA+/f[E/5l] (A1) 


Nous avons 


P°=(p+gA)}+2/f[(E/5).(p+gsA) 
+ 2f2E2+ if (RotE.c) 


ES PROPRIÉTÉS DES CORPUSCULES DE SPIN 1/2 Af?* , 993 


et l'équation (40) devient 


| 


I 


! LEE 3 
HT j SRE 
sh JE? 
es RER RU Fe 
[EE f|Go)+ 2MoC Ya = 0. (42) 


Le Système (42) correspond aux équations d’ondes 
à l’approximation non relativiste cherchées. 
Ces équations introduisent l’impulsion généralisée 


P=p+sA+/[E/o] (43) 
qui correspond à l’approximation non relativiste 
du vecteur — IT et du tenseur d’impulsion 


Lg 
2/M0C 


p+sA+ [HAa]-/[E A ul. 


En effet, ce dernier s’écrit encore 


Kg 
2 MoC 


p+sA— [H À o]p—JflE /\o]e: 
et si nous prenons seulement les termes d'indice 22 
en p, tenant compte de (P2)29 — 0, (05)29 = — 1, 
ce tenseur se réduit à l’expression (43) de P. 

De même, avec 


n=(p+gA)u—f([E À 5]—Ho) 
= (p+sA)p5—/([E As] — Hi) 


les termes d’indices 22 en p nous donnent 


(Ms = (p+£sA+f[E Ac) =—-P. 


L’'équation (42) considérée comme généralisant 
l’équation de Schrôdinger introduit les deux moments 
magnétiques propres en g et f d’une façon additive 
dans le terme d’énergie potentielle 

AO 
| & —/e|(.0) 
l’impulsion généralisée (41) remplaçant l’impul- 
sion p + gA usuelle. 


6. Nous allons examiner plus en détail l’approxi- 
mation de l'équation (6) lorsque le champ extérieur 
se réduit à un champ électrique central E = E (r) 
en mettant en évidence le couplage spin-orbite. 

L’équation (6) se réduit alors à 


| po+ gU(r)+pa—/f{[E(r).zx]+ mo Cas |Y = 0. (44) 
La décomposition (36) nous donne encore en posant 


W 
Po = Po + Mc, Po Ÿ = ee Ÿ, W € moc?, 


e + EURE »moc) o1+ [(p +ifE).5]%:=0, 


W (45) 
(= + £& ü) 29 +[(p —if E).c]y1= 0. 
Posant 
Des RMS Si ET 16 
W + gcU + VRP rl AUS (46) 
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nous en tirons 


c [Cp a ifE).0]%, 


+sU— L (47) 


ue 2 M0 C 
x [(p eo 
Calculant l'expression 


K(p—1f#£).3] (Ip +iÿfE).s] 
=/f()l@—4#ŸE#).31[(p+4/E).01] 


À in (ÆS ) [tp + if E).0] 
le premier terme du second membre s ‘écrit 
fCr)Ip+ fE?+ 2 /[(È À ç).pl 
=fCe+/EAs))=rET] 
=f(r)IP PE) | (48) 
En utilisant les relations vectorielles 
-(A.c)(B.5)=(A:B)—([A/Bl|.c0), 
([x A Alo)=(x.TA À 61), 
(xX.A) =TAr 
le second terme devient 


(x.5)[(p + /E).0] 
=(xp)—({x/pl.is)+1/(&.E)+/([x/AEl.0) 
= (x.(p+/TE À dj) —([x /\ pJ.io)+i/(&.E). 

Or, si E — E (r), nous avons 
(x.[E A sl)=0o 


_et le terme cherché se réduit à 


ie FEI ue. A pls)| 
ce even) (ES a] 


Introduisant le moment orbital 
Lx N\ p} O0) 

l'équation approchée (47) s'écrit maintenant 
Cie er mass 


ce 4 2 MoC 


Prin Rs (Jess «n 


Or, nous avons, en nous bornant aux premiers 
termes, 


P:—/'E:) 


W 


int dr 
> Mo C PATES 2M)C — (2 mc)? 
aa , AO [ œ 21 
4 (OR EE (7 
n(r)= 6 _ _ le] ( ) fe 
2M0C* (2moc} (2moc)? 


ns 
_Introduisant ces expressions dans l'équati 


nous avons les équations approchées aux ie 
ordres É 


D ui Pepe) Le 0 2 


- 2700C 


Nous retrouvons ici l'équation (42). 
20 À l’approximation suivante : © 


1e LE Dr) Re SRE E. SE 
Pal [P2= PE Fr E “2 
gME(r) d }= _ *)] es. 
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Cette expression ne fait intervenir l'interaction 
en / que dans l'impulsion généralisée P, mais non. 
G\ 24 


La 
pas dans le terme de couplage spin-orbite (1 


qui ne dépend que du coefficient g. 4 
Si nous introduisons, au lieu du moment orbital | 


L=[x A pl]. 
le moment orbital généralisé 


= [x sAPI= [RAGE AD]. 


— 


nous avons * 
(£.6) =([x A pl: 9+s([r A 1841] <) Æ Es 
Or, tenant compte de TER 
([E /\ +].x) = 0 

ous avons = 


“([zAIEAS]|.) é 
SEA Æ)—2(x. E)=— (x. E), + 


= RE 2 ou 


pour E=E(7) 


d’où 


Introduisant cette expression dans (62), : cette 
équation devient ; 


bee gU(r)— [P°— fÆ] 


( 2 M0 C 


Crop Pi ER 0 

ghE(r) R.0 | 4 
ne D 1 
mettant en Qu côté du terme de” correction | 


relativiste en ce le terme de couplage Que 
orbite en (#2). ; 
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